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Abstract. In this paper we obtain the asymptotic behavior of solutions of the 
Klein-Gordon equation on Lorentzian manifolds (X°,g) which are de Sitter- 
like at infinity. Such manifolds are Lorentzian analogues of the so-called Rie- 
mannian conformally compact (or asymptotically hyperbolic) spaces. Under 
global assumptions on the (null)bicharacteristic flow, namely that the bound- 
ary of the compactification X is a union of two disjoint manifolds, Y±, and 
each bicharacteristic converges to one of these two manifolds as the parameter 
along the bicharacteristic goes to +00, and to the other manifold as the pa- 
rameter goes to —00, we also define the scattering operator, and show that it 
is a Fourier integral operator associated to the bicharacteristic flow from Kf 
to Y-. 



1. Introduction 

Consider a de Sitter-like pseudo-Riemannian metric g of signature (1, n — 1) on 
an n-dimensional (n > 2) manifold with boundary X, with boundary Y, which 
near Y is of the form 

dx 2 — h 
9 = 2 — ' 

h a smooth symmetric 2-cotensor on X such that with respect to some product 
decomposition of X near Y, X = Y x [0, e) x , h\y is a section of T*Y <g> T*Y 
(rather than merely TyX ® TyX) and is a Riemannian metric on Y. Let the wave 
operator □ be the Laplace-Beltrami operator associated to this metric, and let 
P = P{\) = □ - A be the Klein-Gordon operator, XeR. 

Below we consider solutions of Pu = 0. The bicharacteristics of P over X° are 
the integral curves of the Hamilton vector field of the principal symbol 02 (P) (given 
by the dual metric function) inside the characteristic set of P. As g is conformal to 
dx 2 — h, bicharacteristics of P are reparameterizations of bicharacteristics of dx 2 — h 
(near Y, that is). Since g is complete, this means that the bicharacteristics 7 of 
P have limits lim^-i-oo j(t) in SyX, provided that they approach Y. While many 
of the results below are local in character, it is simpler to state a global result, for 
which we need to assume that 

(Al) Y = Y + U Y_ with Y + and Y_ a union of connected components of Y 
(A2) each bicharacteristic 7 of P converges to Y + as t — > +00 and to Y_ as 
t — * — 00, or vice versa 
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Due to the conformality, the characteristic set S(P) of P can be identified with a 
smooth submanifold of S*X, transversal to dX, so SyX fl S(P) can be identified 
with two copies S±Y of S*Y, one for each sign of the dual variable of x. Under our 
assumptions we thus have a classical scattering map S c \ : S*LY+ — > SZ.Y-. 

It is well-known, cf. 0J, that (Al) and (A2) imply the existence of a global 
compactified 'time' function T, with T G C°°(X), T\y ± — ±1, and the pullback of 
T to S* X having positive/negative derivative along the Hamilton vector field inside 
the characteristic set S(p) depending on whether the corresponding bicharacterstics 
tend to Y + or Y~_ . Notice that 1 — x resp. x — 1 has the desired properties near 
Y + resp. Y-, so the point is that a function like these can be extended to all of 
X. Moreover, such a function gives a fibration T : X [—1,1], hence X is in fact 
diffcomorphic to [—1, 1] x S for a compact manifold S. In particular, Y + and Y_ 
are both diffcomorphic to S. Denote the level set T — to by St . With any choice 
of such a function T, a constant to G (—1, 1), and a vector field V transversal to 
St (e.g. take the vector field corresponding to dT under the metric identification 
of TX° and T*X°), P is strictly hyperbolic, and the Cauchy problem Pu = in 
X°, u\ St0 = -00, Vu[ Sto = ipi, ipo,^i G C°°(S to ) is well posed. 

Theorem 1.1. (See Theorem\6Ji) Let s±(A) = ^ ± v/ ' (n ~ 1)2 - A~. Assuming 
(Al) and (A2), the solution u of the Cauchy problem has the form 

(1.1) w^ar^V+^-^V, w ± eC°°(A), 

«/ s+(A) — S-(A) = 2-y/ ( -"~ 1 - > A is not an integer. If s+(A) — s_(A) is an in- 
teger, i/ie same conclusion holds if we replace w_ G C°°(A) by u_ = C°°(A) + 
x s+{\)-s-(\) \ ogx c co (X). 

Conversely, the asymptotic behavior of v± either at Y + or at y_ can be prescribed 
arbitrarily, see Theorem 15. 51 Thus, assuming Al and A2, if s+(A) — S-(A) is not an 
integer, we show that given g± G C°°(T+) there exists a unique u G C°°(X°) such 
that Pu = and which is of the form (jTTTJ) and such that 

(1-2) v+\y + = 9+, V-\ Y+ =g— 

If s + (A) — s_ (A) is a non-zero integer, the same conclusion holds if we replace V- G 
C°°{X) by w_ =^i=o AKS_(AK1 ai^+a; s +W- s -( A ) \ogxC°°(X), a s G C°°{Y), see 
Theorem 15.51 For A = ^ . ' , a similar results holds, with 

(1.3) u = x^ n - 1 ^ 2 v + + x {n - 1)/2 \ogxv-, v± GC°°(A), v±\ y+ = g±. 

That is, for all A G M, there is a unique solution of Pu = with two pieces of 
'Cauchy data' specified at Y + . Note the contrast with the elliptic asymptotically 
hyperbolic problem (conformally compact Riemannian metrics) : there one specifies 
one of the two pieces of the Cauchy data, but over all of Y (not only at Y+), see 
[9j. The quantum scattering map is the map: 

S : C oc (Y + ) © C oc (Y + ) — » C°°(F_) © C°°(F_), S(g + ,g^) = (v + \ Y _,V-\y_). 

Of course, the labelling of Y+ and T_ can be reversed, so S is invertible. In fact, it 
is useful to renormalize S = 5(A) somewhat so that the two pieces of Cauchy data 
at infinity carry the same 'weight'. Let A'^ denote the operator which is on the 
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orthocomplcment of the nullspace of A^ and is the identity on the nullspace, so A' h 
is positive and invertible. Then the renormalization is 

5(A) 

= ((A' h )- s+{x)/2+n/4 © (A; i )- s -( A ) +n / 4 )5(A)((A , h ) s +( A )/ 2 -"/ 4 © (A' h ) s -( A )/ 2 -"/ 4 ); 

this is analogous to using Aipo in place of ipo for the finite time Cauchy data, where 
A G ^ 1 (S to ) elliptic, invertible. We show that: 

Theorem 1.2. (See Theorem \7.21\ ) Suppose that s + (A) — s_(A) is not an integer, 

i.e. A ^ ( "- 1} 4 - m , in £ N. S = S(X) is an invertible elliptic Oth order Fourier 
integral operator with canonical relation given by S c \, and S is a Fourier integral 
operator. 

Remark 1.3. The somewhat strange powers in the normalization correspond to 
making the map from Cauchy data at infinity to Cauchy data at time to G (—1,1) 
a FIO of order 0; see Proposition 17.201 

Note that the canonical relation is independent of A. While our parametrix 
construction for S(X) does not work apparently if s+(A) — s-(A) is an integer due 
to the possible non-solvability of a model problem with the prescribed ansatz, it 
is expected that with more detailed analysis (changing the ansatz slightly to allow 
logarithmic terms in x) one can prove the theorem in this case as well. Moreover, 
we actually construct a parametrix for the solution operator (g+,<7_) t— > it, and 
even if s + (A) — s~(A) is an integer, the part of the operator corresponding to g+ 
(i.e. with g_ = 0) can be constructed as a Fourier integral operator. 

In addition, if g is even, i.e. there is a boundary defining function x such that 
only even powers of x appear in the Taylor series of g at dX expressed in geodesic 
normal coordinates, see [5] for the Riemannian case, then the log a; terms in v_ 
disappear and our parametrix construction for 5(A) goes through provided that 
s+(A) — s_(A) is odd. In particular, this covers the actual d'Alembertian (A = 0) 
if n is even. 

For the Cauchy problem, we similarly have: 

Theorem 1.4. For t Q G (-1,1) and for all (tp ,ipx) G C 00 ^) 2 , let u G C°°(X°) 
denote the unique solution of the Cauchy problem Pu — in X° , u\g t ~ ip , 
Vu\s tQ — "01- This solution u has asymptotic expansion as in (II. ip . If A ^ 

^ n — — , m £ M, the operators 

(^C^i) »-> (v+\y+,V-\y + ) and (ipo,tpi) i-> (v+ |y_ , w_ \ Y _ ) 

are both Fourier integral operators associated to the bicharacteristic flow. 

To justify our terminology of asymptotically de Sitter spaces, we recall that de 
Sitter space is given by the hyperboloid 2 +1 + 1 in E™ +1 equipped 

with the pull-back of the Lorentzian metric dz^ +1 — dzf — ... — dz 2 . Introducing 
polar coordinates (r, 0) in the first n variables and writing t — z n+ i, the hyperboloid 
can be identified with R t x §g _1 with the Lorentzian metric 

with d0 2 being the standard Riemannian metric on the sphere. For t > 1, say, we 
let x — t _1 , and note that the metric becomes ^ 1+x - — dx ^~' yl+x - d6 , which is of 
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the required form. An analogous formula holds for t < —1, so compactifying the 
real line as an interval [— 1, l] s (with s = 1 — xforx<^, say), we see that de Sitter 
space indeed fits into our framework. (Thus, one can take T — s for the global 
compactified time function.) We also note that another, perhaps more familiar, 
form of the metric can be obtained by letting t = sinh p; the metric becomes 
dp 2 — cosh 2 p d6 2 . (One can take e.g. T = tanh p here.) 




Figure 1 . On the left, the compactification of de Sitter space with 
the backward light cone from q + and forward light cone from 
are shown. f2 + , resp. f2_, denotes the intersection of these light 
cones with T > 0, resp. T < 0. On the right, the blow up of de 
Sitter space at g+ is shown. The interior of the light cone inside 
the front face S q+ can be identified with the spatial part of the 
static model of de Sitter space. 



We also use this occasion to explain the connection with the static model of 
de Sitter space. This corresponds to singling out a point on Sg -1 , e.g. qa = 
(1,0,..., 0) e C K". The static model of de Sitter space then is the in- 

tersection of the backward lightcone from qq considered as a point q + on Y + (so 
T(q + ) = 1) and the forward light cone from qo considered as a point g_ on Y_ (so 
T(g_) = —1). These happen to intersect the equator T = (here t = 0) in the 
same set, and altogether form a 'diamond', see Figure [1] Explicitly this region is 
given by z\ + . . . + z\ < 1 inside the hyperboloid. The standard static coordinates 
(r, r, u>) on the 'diamond' are given by 



- ■ l z% + ... + zl = x l + z 2 



n+l 



sinliT = 



Zn+l 



1 *n+l 
-li 



LU = r- 1 (z 2l ...,z n )eS n - 2 . 

In these coordinates the metric becomes (1 — r 2 ) dr 2 — (1 — r 2 )~ x dr 2 — r 2 dui 2 . Note 
that the singularity at r — 1 is completely artificial (is due to the coordinates), the 
metric is incomplete, but is conformal to a complete Lorentzian metric, of product 
type, with □ also of product type. While one can analyze the solutions of the 
wave equations on de Sitter space at points inside the 'diamond' by considering 
the diamond only (in view of the finite propagation speed for the wave equation), 



THE WAVE EQUATION ON ASYMPTOTICALLY DE SITTER-LIKE SPACES 5 

the resulting picture does include rather artificial limitations. For instance, the 
asymptotics at the sides of the diamond are automatically smooth in de Sitter 
space (as we have a standard wave equation there), which is not obvious if one's 
world consists of the diamond, and the local static asymptotics, corresponding to 
the tip of the diamond at Y+, describes only a small part of the asymptotics of 
solutions of the Cauchy problem on de Sitter space. However, the 'spatial' part of 
the static operator (or modifications of it) do show up in our analysis as models 
for the Poisson operator (<? + ,<?_) i— ► u; the proper place for its existence is on the 
interior of the light cone in the blow up of the distinguished point q+ in de Sitter 
space. 

It should be pointed out that the de Sitter-Schwarzschild metric in fact has many 
similar features, and the analogous result is the subject of an ongoing project with 
Antonio Sa Barreto and Richard Melrose. Weaker results on the asymptotics in 
that case are contained in the part of works of Dafermos and Rodnianski concerned 
with the underlying linear problem [2] (they study a non-linear problem) , and local 
energy decay was studied by Bony and Hafner 1 , in part based on the stationary 
resonance analysis of Sa Barreto and Zworski |16j . 

We also note that on de Sitter space itself, one can solve the wave equation 
explicitly, see [T5], but even the 'smooth asymptotics' result, Theorem 16.11 is not 
apparent from such a solution. 

There are two rather different techniques used to prove the results here. The 
'rough' results yielding the existence of the asymptotics, Theorems 15. 51 and 16. 1[ are 
proved using positive commutator estimates, which roughly speaking describe the 
microlocal (i.e. phase space) propagation of L 2 (or Sobolev) mass ('energy'). Such 
methods are very robust, but (unless they are used in a more sophisticated form as 
in [6]) give less precise results. The Fourier integral operator results are proved by 
a parametrix construction which is significantly more delicate (taking up two-fifth 
of this paper), but is very instructive. It is at this stage that the static de Sitter 
model shows up on the front face of [X x Y + ; diag y+ ] ; see P a in Section [7] One 
should think of this as analogous to the way the hyperbolic Laplacian shows up as a 
model on the front face of the 0-double space for conformally compact Riemannian 
manifolds, see [9]. 

The plan of the paper is the following. In Section [2] we adopt a O-microlocal 
point of view, and analyze propagation of singularities in the O-cotangent bundle 
introduced by Mazzeo and Melrose [9] two decades ago. The proof uses positive 
commutator estimates, and is quite similar to propagation of singularities for man- 
ifolds with boundary equipped with a so-called (incomplete) edge metric, which 
includes e.g. manifolds with conic points - see [14] and [13] and references therein. 
In the following sections we analyze local solvability near the boundary as well as 
conormal regularity of the solutions there. We emphasize that the results of Sec- 
tions'^ do not need the global assumptions (A1)-(A2). In Section [5] we prove a 
unique continuation theorem at dX (i.e. at 'infinity') by a Carleman-type estimate, 
and use it to prove that the asymptotic behavior of the solutions in fact determines 
the solutions, i.e. we can talk about a 'Cauchy problem at infinity', hence also 
about the scattering map. In the final section we construct a parametrix for the 
scattering map, and use it to show that it is indeed a Fourier integral operator. 

I am very grateful for Rafe Mazzeo, Richard Melrose, Antonio Sa Barreto and 
Maciej Zworski for numerous fruitful discussions. In particular, I thank Richard 
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Melrose for pointing out that the assumptions (Al) and (A2) imply the existence 
of a global time foliation, while relating the analysis here to the static de Sitter 
model arose from discussions with Maciej Zworski. 

2. 0-GEOMETRY AND PROPAGATION OF O-SINGULARITIES 

For the purposes of analysis, we need a good understanding of bicharacteris- 
tic geometry. Thus, note that P G Diff (X), in the zero-calculus of Mazzeo and 
Melrose [9]. Let °T*X denote the zero-cotangent bundle of X. Its elements are 
covectors of the form £ — + rj Then the principal symbol p = cr(P) is a homo- 
geneous degree 2 polynomial on °T*X; explicitly at Y, p\y — £ 2 — H\y, H\y the 
metric function corresponding to h, and p itself is the metric function of the dual 
pseudo-Riemannian metric g. We refer to [9l [16] for nice descriptions of the basic 
setup, and [T4J, [T3] for analysis of a hyperbolic problem in the related edge setting. 

If a is a homogeneous function on °T*X \ o, then there is a (homogeneous) 
Hamilton vector field H a associated to it on T*X° \ o. A change of coordinates 
calculation shows that in the O-canonical coordinates given above 

H a = (<9 5 a)(£<% + r)d n + xd x ) + x(d v ad v - d y ad v ) - ((xd x + + rjd v )a)d 5 , 

so H a in fact extends to a C°° vector field on °T*X \ o which is tangent to °Tg X X. 
At x — this gives H a = {d^a)R — (Ra)d^, where R is the radial vector field 
£9^ + 7]d v on °T*X. Since a is homogeneous of degree, say, k, Ra = ka, and 
H a = (d^)aR — kad^, so on the characteristic set S(a) = a _1 ({0}) of a, at x = 0, 
H a is radial. It is thus rather convenient to consider the cosphere bundle °S*X 
which is the boundary at fiber infinity of the fiber radial compactification a T*X of 
^X 1 * X 

As we work with p, so that near Y, £ ^ on the characteristic set, we use 
projective coordinates rj = r)/\£\, p = valid near S(p). Then 

(sign £,)p^ 1 H a = - ((pd p + f)df))a) {-pd p + xd x ) + x(df,ad y - d y adf,) 

+ ((xd x - pd p )a) {pd p + , qdf l ). 

Thus, for a G p- k C° c ( f*X), W a = p^Ha is a smooth vector field on °f*X, 
whose restriction to °S Y X is (sign ^kafidfj, i.e. it vanishes at a = 0. Thus, if da 
is not conormal to °S' Y X in °S*X, so S(a) is transversal to °S Y X, then W a is 
a smooth vector field on E(a) that vanishes at x = 0, and hence is of the form 
W a = xW' a , W' a G V(S(a)). 

Applying this with a = p = p- 2 C°°( T*X) yields that inside S(p), W p = xW' p , 

^;U=o= (sign £){2d x + H h ), 

Hh the Hamilton vector field of h (evaluated at (y,?))). In particular, Wp is 
transversal to Y. Also, W p is complete, and 7 is an integral curve of W p , then 
a reparameterized version of 7 is an integral curve of W p , hence lim t ^_( s ; gn ^) 00 j(t) 
exists in °S Y X n S(p). (Note that the map switching the sign of covectors pre- 
serves even functions, such as p, while transforms H p to —H p .) Conversely, for 
any q G °S Y X n S(p) there is a unique (up to translation of the parameteriza- 
tion) integral curve of W p with limit q as t — > — (sign£)oo, namely this is just a 
reparameterization of the unique integral curve of W p through q. Note also that 
°S Y X n S(p) can be identified with two copies of S*Y, one for each sign of £; we 
write these as S+Y and S1Y. 
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Suppose now that Y = Y+ U Y— , where Y± are unions of connected components 
of Y, and this decomposition satisfies that all bicharacteristics 1 1— > 7(t) of P satisfy 
lim t ^ +oc 7(t) e 5*1+, limt_»_oo 7(t) G S'*F_, or vice versa, i.e. that (Al) and (A2) 
of the introduction hold. For q G SI. 1+ there is a unique bicharacteristic with 
lim t ^_ oc 7(t) = g. By (Al) and (A2), \im t ^ +00 ^f(t) = q' G S'*l r _ exists; as we 
saw above, it necessarily lies in S!¥—. The classical scattering map is the map 
S c i : S*Y + — > with <S c i : q ^> q' . Fixing a homogeneous degree 1 function on 

T*Y \ o, we can extend these to maps T*Y + \ o — > T*Y"_ \ o - we can use /i 1 / 2 , for 
instance. The induced relation on (T*Y + \o) x (T*Y"_ \o) is Lagrangian with respect 
to the twisted symplectic form (i.e. with a negative sign on one of the factors). 

As follows easily from the results of [4], (Al) and (A2) imply the existence of 
a global compactified 'time' function T, with T G C°°(X), T\y ± = ±1, and the 
pullback ir*T of T to S* X having positive/negative derivative along the Hamilton 
vector field inside the characteristic set E(p) depending on whether the correspond- 
ing bicharacterstics tend to Y + or Y_ . Notice that 1 — x resp. x — 1 has the desired 
properties near Y"+ resp. Y_ , so the point is the interior of X these can be extended 
to all of X. 

With any choice of such a function T, a constant to G (— 1, 1), and a vector field 
V transversal to St P is strictly hyperbolic, and the Cauchy problem Pu = in 
X°, u\ StQ = -00, Vu\ Stg = Vi, Vo, ipi G C oo (S' t0 ) is well posed. 

Our first result is that of O-regularity of solutions of Pu = with a weight given 
by a space u a priori lies in. There is a dichotomy between solutions depending on 
the a priori regularity relative to this weighted space. If the a priori regularity is 
low, we only obtain regularity up to a limit implied by the weight, but we do so 
without having to assume any interior regularity for u. If the a priori regularity 
is high, then we obtain additional regularity up to the limit corresponding to the 
smoothness of u in X°. 

Proposition 2.1. Suppose that q 6 7, and suppose that u is in Hq°' s °(X) in a 
neighborhood of q and Pu = 0. Then: 

(i) If 7'o < Sq + 1/2 then u is in Hq' s ° (X) near q for all r < s + 1/2. 

(ii) Ifr > s + 1/2 andr> r , a G °S*X n S(p) then a WF r Q ' s ° (u) provided 
that the bicharacteristic 7 approaching a is disjoint from WF r (u) C S*X° . 
The same conclusion holds if ro < sq + 1/2, but a ^ WFj 1,s °(«) for some 
n > so + 1/2. 

(iii) In particular, if ro > So + 1/2 and r > ro, then u is in Hq S °(X) near 
q provided that all bicharacteristics approaching °S*X are disjoint from 
WF r (u) C S*X°. 

Proof. This proof is very similar to the proofs of propagation of 'edge regularity' 
for the wave equation with incomplete metrics in [14] and |13j , so we shall be brief. 

While pH p restricts to a smooth vector field on E(f>) with vanishing restriction at 
Y, if we evaluate pH p as a section of the b-tangent bundle of °T*X at E(p) n°SyA, 
the result is more interesting: pH p = 2(—pd p + xd x ) in this sense. Correspondingly, 
for A G * n ''(A), the symbol of i[P,A] G ^ +1 ' l (X) is H p a = 2(m + 
a = cr(A), at S(p) n SyX. Thus, much as [M] and [15] . one can show propagation 
of zero- regularity into the boundary for m + Z 7^ 0. Unlike in the setting of [14] . 
the characteristic set of P only intersects the boundary Y in radial points, i.e. 
there is no propagation inside Y , which explains why there is no requirement for 
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m + I having a particular sign (as long as it is non-zero), although the results are 
different depending on the sign: (i) has no wave front set assumptions on u. This 
corresponds to the presence of a cutoff x, identically 1 near Y, such that d x \ < 0, 
the sign of the commutator with \ agrees with the sign arising from the weights 
if m + / < 0. Moreover, one can microlocalize in S Y X by pulling back functions 
from SyX (~l E(p) using the flow of Wp, extending them to a neighborhood of the 
characteristic set in an arbitrary smooth fashion. 

Thus, let G C°°{S Y X n E(p)), and for any integral curve 7 of W p with 
7(0) e^nS(p)),we let V(7(*)) = ^o(7(0)). Note that this defines a C°° function 
on E(p) near F, for the map $ : 5f A" n E(p) x [0, e) -> E(p) given by $(t?,t) = 
exp(tVTp)(7 is a local diffeomorphism near i = 0. As E(p) is a C°° submanifold of 
S'*A, we can extend ip to a C°° function on still denoted by ip, hence further 

to an element of C°°(°f*X), at least near Y. Now let \ G C~([0, e)) be such that 
x' = — X07 X = 1 near 0, x > 0, x 1 / 2 is C°°, and let a — p~ m x l x(x)ip, and note that 
W P x(x) = b 2 xx'{x) with b > near F. As vanishes at p = 0, we deduce that 

ff p a = 2(m + l)p- m -Vx(!B)V' + 2p- m - V+^YM + pp~ m+1 x l e + p~"V/, 

with 6, e, / G C°°(°T*A). 

Now the standard positive commutator argument finishes the proof of the propo- 
sition, see e.g. [14] . For the reader's convenience, we sketch the argument, skipping 
the (necessary but straightforward) regularization part of the argument. Thus, 
c-to-i(*[-Pj A]) = H p a shows that 

i[P, A] = 2(m + l)AA*AA - B*B + PE + F, 
A6$f' V2 (I), a(i) = a(A) 1 / 2 , 

B G ^ m+1)/2 '' /2 (A), WF'(B) C suppxo nsuppV', <r{B) = b X o(2p- m - V +1 ) 1/2 
E G *o n_1,/ (A), F G *o M (A), 

A G ^0 2 '°{X) elliptic formally self-adjoint with positive principal symbol, p^ 1 / 2 . 
Proceeding as in [14] shows that for u with Pu = 0, 

I ± \\AAu\\ 2 H o,o (x) - \\Bu\\ 2 H o. 0{x) \ < C\\u\\ 2 K/2 ^ l/2(xy 

provided that the right hand side is finite, with the — sign applying if m + I < 0, 
and the + sign applying if m + I > 0. In the first case, the second term on the 
left hand side can simply be dropped, so we do not need to make any assumptions 
on the ij( m+1 )/ 2 norm of u, while in the second case we need to assume that 
WF' m+1 " 2 (u) is disjoint from suppxo, m order to conclude that ||AA.ti|| fl -o,a^ 

is finite, i.e. WFQ m+1 ^ 2 ' ^ 2 (u) is disjoint from the elliptic set of A, i.e. from the 
interior of suppi/; near x = 0. 

The standard iteration argument now proves the proposition. □ 

The approximation process prevents us from crossing the line r = so + 1/2, which 
is why we cannot directly obtain information about u in Hq S ° (X) with r > sq + 1/2 
unless we know u is in Hq°' s ° (X) for ro > sq + 1/2. However, if u G H r °' s ° (X) with 
r = s + l/2-e/2, so r < s + l/2, then u G ff r °' s °- £ (A), and r > (s -e) + l/2 
now. We thus deduce: 
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Corollary 2.2. Suppose that q £ Y , and suppose that u is in Hq°' s °(X) in a 
neighborhood of q and Pu = 0. If r > ro and s < so, a £ °S*X n S(p) then 
a ^ WFq S (u) provided that the bicharacteristic 7 approaching a is disjoint from 
WF(m) c S**A°. 

in particular, u is in Hq S (X) near q provided that all bicharacteristics approach- 
ing °S*X are disjoint from WF r (u) C S*X°. 

Remark 2.3. Thus, we gain full 0-regularity for solutions if we are willing to give 
up some (arbitrarily little) decay. Note that (ii) of the Proposition states that one 
can take s = sq if r<j > sq + 1/2, so the present corollary is only interesting if 
r < s + 1/2. 

Proof. Let s < sq be given, and let e = sq — s > 0. As remarked, we may assume 
ro < so + 1/2, and if needed, we can decrease r so that r < s + 1/2. By (i) of 
Proposition O a <£ WF r ' s °(u) for all r < s + 1/2. Then a <£ WF s 0+1/2 ~ e/2 ' s ° (u), 
and hence a <£ WFo° +1/2 " e/2 ' s °" e (u). By (ii) of PropositionED a <£ WFj s »" f («) = 
WFq s (u) for all r, proving the corollary. □ 



3. Local solvability near dX 

In this section we show the solvability of Pu = near dX in suitable senses, 
P = □ — A. This relies on a positive commutator estimate with compact error term, 
so we need to control the normal operator of our commutator in the 0-calculus. 
Recall from [5] that the normal operator map on Diff^X) (or ^q(X)) captures 
Q £ DxSq(X) modulo xDxSq(X), as opposed to the principal symbol map, which 
captures it modulo Diff -1 (X). The compactness referred to above then is that of 
the inclusion map for the associated Sobolev spaces, Hq S (X) to Hq ' s (X), with 
r > r' , s > s'l note that compactness requires improvements in both the regularity 
and decay orders, hence control of both the principal symbols (described in the 
previous section) and normal operators. 

We thus start by calculating the normal operator of P, as well as that of its 
commutator with another operator A. Thus, we calculate the the commutator 
modulo terms with an additional order of vanishing. As P £ Diff^A), and our 
commutant will be an operator A r £ x r ~ x Diff (A), [P, A r ] £ x r_1 Diffj^A), so we 
need to compute [P, A r ] modulo x r Diff^A). This is computation is thus unaffected 
if P is changed by addition of a term in a;Diff (A), or A r is changed by a term 
in x Diffg(A). This means that effectively we may assume that A has a product 
decomposition near Y and h is actually a Riemannian metric on Y. The wave 
operator is the Laplace-Bcltrami operator associated to this metric: 

□ = (xD x ) 2 + i(n - l)(xD x ) - x 2 A Y = (xD x )*(xD x ) - x 2 A Y , 

with the adjoint taken with respect to the pseudo-Riemannian density x~ n \dxdy\. 

We remark here that the actual normal operator in the 0-calculus (which results 
from restricting the Schwartz kernels to the 0-front face) is even simpler than this 
model, for it localizes in Y. Thus, one could simply compute with the Euclidean 
Laplacian in Y, but as this has absolutely no impact on our considerations, we use 
our more global model. 
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We let A r = x r D x + x r 1 , which is symmetric, and compute 

[P, A] = [{xD x f + i(n - l){xD x ), x r D x + x^ 1 ] 
- [x 2 ,x r D x + ar r_1 ]A y 
= -2t |(r - l){xD x + i T ^-y x r -\xD x + i^-) + x r+1 A Y } . 

Thus, up to the factor —2i, this is clearly a positive operator for r > 1. We would 
like to improve this statement, and in particular show that this is greater than 
Cx r ~ x for suitable C, at least in a range of r, and at least modulo terms of the 
form PB + B*P. 

The flexibility we have here in arranging this positivity is the choice of the 
coefficient B of P. Thus, we convert part of the tangential Laplacian term, x r+1 Ay 
into P by writing x r+1 Ay = fx r+1 Ay + (1 — 7)i r+1 Ay, with 7 to be determined, 
and writing 

x r+1 A Y = ^{x r -\{xD x )\xD x ) -\-P) + ((xD x )*(xD x ) - A - P)x r ~ 1 } 
in the first term. We deduce with B = — 3r x r ~ 1 , 

l -[P, A] =(r - l)(x^ + iV^Ly x '-\ x D x + i^I) + (1 - 7 )z r+1 A y 

+ ^x r -\xD x )*{xD x ) + l^)'^)^- 1 - 7 Aa; r - 1 + P-B + P*P. 

Now, the form of the first term is quite convenient to us in view of the factor x r ~ 1 , 
corresponding to a weighted estimate on a; _ ( r_1 )/ 2 -L 2 relative to x~ n dx, since its 
null-space consists of x( n ~ r ^ 2 , which just misses being in x^^^ 1 ^ 2 L 2 (i.e. is in 
x -( r - 1 )/ 2 - s L 2 f or a u g > g0 ^ w jn gj ve us optimal zeroth order terms below, 
and saves us having to use that for all s, 

(3.1) {2s - n 4 - 1)2 \\x^u\\ 2 <\\^D x u\\ 2 . 

Note, however, that the first term can easily be written in a simpler looking form, 
(xD x + i 1 ^L)*x r - 1 {xD x + = {xD x yx r -\xD x ) - ( n ~ r ) ^-1 

This can be checked easily as the two sides have the same principal symbol, so 
their difference is first order, moreover both sides are real and self-adjoint, hence 
actually zeroth order, i.e. multiplication by a smooth function. Their equality can 
be checked by evaluating them on 1. Moreover, a similar calculation yields 

^(x^ixD^ixD^ + (xD^ixD^- 1 ) 

1 .n-r. t .n-r. (n + r-2)(n-r) __, 
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Thus, 



(3.2) 



- [P, A] =(r - 1 + 1 ){xD x + i-—)* x r-\ xDlB + i— )) + (1 - 7 K' +1 < 



(n — r)(n + r 



+ 7 ^— " " 4 ' " 1 - AJ x 1 - 1 + PS + P*P. 

In order to obtain a 'positive commutator', modulo the terms involving P, we thus 
need that 

/ (n — r)(n + r — 2) 
(3.3) r-l + 7, 1-7 and 7 ( — A 

* in— r)(n+r — 2) \ n 

As 1 i - A = gives 





we introduce 

(3.4) 1(A) = Re „ , 

so /(A) = for A > Z(A) > for A < 

Lemma 3.1. The quantities listed in (|3.3[) have the same (non-zero) sign if: 

• if r > max(0, 1 — 21(A)), r =/= 1 + 2Z(A), in which case they are all positive, 
or 

• if r < min(0, 1 — 2Z(A)), in which case they are all negative. 

Proof. First, note that for G (-l(A), 1(A)), i.e. r G (1 - 21(A), 1 + 2/ (A)), 

(n-r)(n+r-2) _ ^ > ^ for ^-1 g [.J^y^ (n-r)(n+r-2) - A < 0. 

For r > 1, r ^ 1 + 21(A) it is easy to arrange that all three quantities in (|3.3p 
have the same sign since the first two terms are positive if I7I is sufficiently small, 
so choosing the sign of 7 correctly, the last term can also be made positive as long 
as r ^ 1 + 21(A) (r > 1 rules out r = 1 - 21(A)). 

In general, the first two terms have the same sign if 7 G (1, 1 — r), resp. 7 G 
(1 — r, 1), depending on whether r < 0, resp. r > 0, and this sign is negative, resp. 
positive in the two cases. 

Suppose first that A < &-jp-. 

If r < 0, we have 7 > 1 by the previous remark, so we need (n + r — 2)(n — 
r) — A < 0, i.e. r ^ [1 — 21 (X), 1 + 2/(A)], which in view of r < amounts to 
r < 1 — 21(A) (and r < 0). In the latter case, if r G (0, 1], 7 > still, but now we 
need (n + r - 2)(rc - r) - A > 0, i.e. r G (1 — 2/(A), 1 + 2/(A)). As r G (0, 1], this 
means r G (max(0, 1 — 2/ (A)), 1]. On the other hand, if r > 1, we have already seen 
that 7 ^ ("-r)(Tt+r-2) _ -A can k e mac [ e positive as well as long as r ^ 1 + 21(A). 

( — 1) 2 

This completes the proof of the lemma if A < {n 4 ' . 

For A > i 1 ^-) 2 , ( "~ r)( " +7 " 2) - A < for all values of r. The 'positive' com- 
mutator criterion thus becomes that r — 1 + 7, 1 — 7 and —7 must have the same 
sign. The first two give 7 G (1, 1 — r), resp. 7 G (1 — r, 1) depending on r < or 
r > 0, as beforehand, while the last two give 7 ^ [0, 1]. As (1, 1 — r) or (1 — r, 1) 
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intersects the complement of [0, 1] in a non-empty set if r < or r > 1, we get ex- 
actly the range stated in the lemma, taking into account that max(0, 1 — 2Z(A)) = 1, 
min(0, 1 - 2Z(A)) = 0. □ 

If the conditions of Lemma [3~T1 are satisfied, the right hand side of (|3.2p . applied 
to v supported near Y, is, modulo the terms involving P, bounded below a positive 
multiple (if all quantities in ()3.3j) are positive), resp. bounded above by a negative 
multiple (if all quantities in (|3.3|) are negative), of the squared x 1 Hq norm of v, 
I = — r rp-. We thus have: 

Lemma 3.2. Suppose 

(3.5) I e (-oo,min(-,Z(A)), I ^ -f(A) or I e (max(-, l(X)), +oo). 

Then there exists C > and 5 > such that 

(3.6) \\x- l v\\Hl<C\\x- l Pv\\ L *. 
for all v G C°°(X) with suppw C {x < S}. 

Remark 3.3. Note that (near x = 0) x s G x'i 2 if / < s — (n — l)/2, so (neglecting 
the \ above) the two critical values I — —l(X) and I — l(X) arise from the monomials 
2;-K A )H — 5 reS p. x'^^+t - , which are exactly the monomial solutions of Pv = 0. 

Proof. Note that (|3.5|) holds if and only if one of the conditions in Lemma I3TT1 holds 

First, suppose that v G C°°(X) supported in x < 5 and g is an exact warped 
product Lorentzian metric for x < 26. Then 

C-Av,Pv) - C-Pv,Av) = C-[P,A]v,v) 

= (r - 1 + ^Wx^ixDx + *^>|| 2 + (1 - ^Wx^dyvf 

(n — r)(n + r — 2) \ r-i „ , , 
^ - A J \\x— v\\ 2 + (Bv,Pv) + (Pv,Bv), 

so as the three squares on the right hand side have coefficients with the same sign, 
\\x- l v\\ 2 H i < C\\x- l Pv\\ L 2(\\x l Av\\ L 2 + \\x l Bv\\ L2 ) 

< Ce^Wx^Pvf^ + Ce(\\x l Av\\ 2 L2 + ||a; ( Pi;||| 2 ). 

As \\x l Av\\ 2 L2 + \\x l Bv\\ 2 L2 < C'Hx-'uH^!, for B = x~ 21 , A = X - 2l {xD x + i^), 
for e > small we deduce that (with a new C > 0) 

\\x- l v\\ H i <C\\x- l Pv\\ L *. 

This proves the lemma for warped product g (with S > arbitrary, as long as on 
x < 25 the metric is warped product). 

If we do not consider an exact warped product metric near Y, then P = Pq + Pl, 
Pq = Do is the wave operator for the warped product metric and Pi G xDiSq(X). 
Moreover, making A self-adjoint with respect to the new metric, A — Aq + Ai, 
A x G a; r DiffJ(A). Thus, 

[P, A] = [P ,A ] + R', R' G x r Di$ 2 (X). 
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Taking into account that I = — i -jr L j for functions v supported in x < S this gives 

\(v,R>v)\<C6\\x- l v\\ 2 H i 

with C depending on R! only (i.e. independent of 6 G (0, 1]), so for sufficiently small 
S > 0, GC1 still holds. □ 



The estimate (|3.6[) gives, by duality, an existence result. As the argument is 
local near each connected component of Y , we have: 

Proposition 3.4. Suppose g is asymptotically de Sitter like, P = □ — A ; l(X) is 
given by (|3.4p , and 

(3.7) I G (-oo,-max(i,/(A))), or I G (- min(-, /(A)), +oo), Z ^ /(A). 

for every f G x l L 2 (X) there exists u G a; Hq(X) such that Pu = f near Y. 
Moreover, if Yj is a connected component of Y , and supp / is disjoint from other 
components of Y , then supp it may be taken disjoint from other components ofY. 

Proof. Note that P = P* (formal adjoint). The result is standard then, see [H 
Proof of Theorem 26.1.7]. Indeed, shows that for / G x~ l H^, v G C°°(X) 

supported in x < 8, 

\{f,v)\<C\\x- l Pv\\ L *. 
Thus, Pv i— > (f,v) is an anti-linear functional on elements of C°°(X) supported in 
x < S, continuous with respect to the cc'i 2 -norm. By the Hahn-Banach theorem 
it can be extended to a continuous conjugate- linear functional on x l L 2 , so there 
exists u G x~ l L 2 such that (f,v) = (u,Pv), and u is now the desired solution for I 
as above. □ 

In order to use the positive commutator argument with v not supported near 
Y, we need a cutoff x> so instead of A — A r , we would really use A — x(x) 2 A r + 
A r x{x) 2 , x = 1 near 0, x £ C£°(R). We can also localize at any given connected 
component of Y; as this can be done by a locally constant function on suppx, we 
do not indicate this in the notation as it leaves the commutator unchanged. Then 

i n — v Ti — v 

- [P, A] =(r - 1 + i){xD x + l —y x r-^(xD x + z^-) 
< 3 - 8 ) + (1 - 7K + VA, + 7 f ^-r)(n + r-2) \ ^ 

+ (xD x )*(x 2 )'(xD x ) + R + PB + B*P, 

where R = R(x), R G C^°(R), supported away from 0. (Again, this comes from 
a principal symbol computation, which has to be carried out away from dX, and 
reality plus self-adjointness shows that R is 0th order.) Thus, modulo the 0th order 
term supported in the interior and terms involving P we have a global 'positive 
commutator' estimate (all terms have the same sign) if r < min(0, 1 — 2/ (A)); if 
r > max(0, 1 — 22(A)) but r ^ 1 + 2Z(A), the commutator terms with y 2 has opposite 
sign compared to the 'main' terms. 

One can also add a regularizing factor, y^f^j = (1 + ex -1 ) -5 with s > small. 

For e > 0, this is a symbol of order — s (i.e. decaying as x — > 0), and is uniformly 
bounded as a symbol of order 0. Moreover, 

(xd x ) k (l + ex- 1 )- 8 = s(l + ex' 1 )'* f k , e , s , 
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where fk,e,s is a symbol of order 0, and is uniformly bounded as such a symbol. 
Consequently, as long as one has a positive normal operator for the commutator 
of P with some operator A, one will also have a positive normal operator for the 
commutator of P with (1 + ex~ r )~ s A(l + err -1 ) -15 if s is small. It is actually 
even easier to simply apply our previous estimate, (|3.6p . to a regularized version 
v e = (1 + ex~ 1 )~ s v of v, for Pv t = (l + ex~ l )~ s Pv + [P, (1 + ea; _1 ) _s ]D, noting that 
(1 + ea; _1 ) s [_P, (1 + ez -1 )^] is bounded by C s in Diff c (X) (c denotes conormal 
coefficients, but should be changed), so the L 2 norm of [P, (1 + ex~ l )~ s ]v can be 
absorbed into the left-hand side of (|3.6[) for s > small. Applying this iteratively, 
we deduce the following: 

Proposition 3.5. Suppose g is asymptotically de Sitter like, P = □ — A, A e K. 
Suppose that u G x 1 "Hq(X) Pu G x l L 2 (X), I > Iq. Suppose also that one of the 
following conditions holds: 

(i) I < -/(A), 

(ii) l > max(±,Z(A)), 

(iii) l > -l(X), I < min(i,/(A)). 

Then u G x l H^(X). 

Moreover, the result is local near each connected component ofY. 

This immediately gives that if a solution of Pu = decays faster than a border- 
line rate, given by ' L 2 , then it is Schwartz. In fact, later in Proposition [5T3J we 
show that such u is necessarily identically 0. 

Corollary 3.6. Suppose that u G x l H${X), k G M, A G M, I > max(i, l(X)), 
PueC°°{X). ThenueC°°{X). 

If the assumptions hold near a connected component of Y only, so does the 
conclusion. 

Remark 3.7. The assumption / > max(i,/(A)) is probably not optimal if 1(A) < ^, 
cf. Remark 15.21 one expects I > l(X) simply. However, this makes no difference 
in the present paper. Moreover, for □ itself this is not a restriction as n > 2 so 
> \. 

This corollary also states in particular that for / G C°°(X) the solution u G 
x 1 Hq(X) of Pu — f near Y, whose existence is guaranteed by Proposition 13.41 is 
in fact in C°°(X). 

Proof. First, we may assume k = 1. Indeed, if k < 1, then I > 1/2 gives k < 1 < 
I + 1/2, so (i) of Proposition 12 . 1 1 applies and gives u G Hq' 1 (X). 

By Proposition 13 . 51 u G x Hl(X) for all I. Thus, by Proposition [HU part 
(i), u G Hq S (X) for all r and s with r < s + 1/2, hence for all (r,s). (Given 
(r, s), consider (r, s') with s' > max(s,r — 1/2) to see that u G Hq S (X) hence 
u G Hq S (X).) In particular, x m Qu G L 2 (X) for all m and all Q G Diff(X), proving 
the corollary. □ 

4. Conormal regularity 

While Proposition 13.41 gives the correct critical rates of growth or decay for 
solutions of Pu = 0, and Corollary 12.21 gives their optimal smoothness in the 0- 
sense, this is not optimal: solutions of Pu = which are C°° in X° are conormal 
to the boundary, i.e. stable (in terms of weighted L 2 -spaces) under the application 
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of b-differential operators. In fact, as usual, cf. [18] and [13], it is convenient to 
work relative to O-Sobolev spaces, i.e. to work with Diffg '5™ (A). However, rather 
than using positive commutator estimates as in these papers, we rely on an 'exact' 
commutator argument (exact at the level of normal operators), much like in |12[ 
Section 12]. Although it was not discussed explicitly in [12] for reasons of brevity, 
the analogous space of operators in that setting would be Diffg C \P C (X), with ty c (X) 
standing for cusp pseudodifferential operators. (Instead, in (12j 'tangential elliptic 
regularity' was used.) 

Definition 4.1. Elements of Diffg Vg(X) are finite sums of terms QA, Q G 
Diffg(X), A G *™(X). We also let x r Diffg be the space of operators 

of the form x r B, B G Diffg V™(X). 

Remark 4.2. Directly from the definition, Diffg &™(X) is a C°°(A)-bimodule (un- 
der left and right multiplication), so in particular x r Diffg i£™(X) is well-defined 
independent of the choice of a boundary defining function x. 

The key lemma is: 

Lemma 4.3. For Q G Diffg (X), A G ^T(X), there exist Qj G Diffg (A), A j G 
&™(X), j = 1, . . . , I, such that QA — A jQj- (With a similar conclusion holding, 
with different Aj, Qj, for AQ.) 

Proof. It suffices to prove the statement for Q G Vo(A); the general case then 
follows by an inductive argument. As V (X) C V b (A), [Q,A] G V™(X), so QA = 
AQ + [Q,A] gives the desired result. □ 

Corollary 4.4. Diffo^b(A) is closed under composition: if A G Diff i&™(X) and 
B G Diffg' ®f(X) then AB G Diffg +fe ' *™ +m '(A). 

We also need the corresponding result about commutators. 

Lemma 4.5. Moreover, if A £ x r ^(X), Q G Diffg (A) then 

[Qi A] G x r Diffg -1 &™(X). 

If in addition a b , m (A)\ bT , 9x = then [Q,A] G x r Diffg *™ _1 (A). 

Remark 4.6. h T*dX is a well-defined subbundle of Tg X X. If we write b-covectors 
as a ^ + r\ ■ dy, then h T*dX is given by x = 0, a = in b T*A. 

Proof. Again, it suffices to prove the first statement for Q G Vo(X). As Vo(A) C 
Vb(X), [Q,A] G "^(A), giving the result for such Q. Iterating this also proves 
that for Q G Diffg (A), [Q,A] G Diffg" 1 V£(X). 

To have the better conclusion, it again suffices to consider Q G Vo(A). As above, 
[Q,A] G *™(X). But, with a = a b , m {A), q = <7 M (Q), 

i(r b<m ([A, Q}) = H a q 

= (d a a)(xd x q) - (xd x a)(d a q) (i d rij a ){ d Vil) _ { d Vi a ){ d Vi'l)) ■ 

This vanishes at h T*dX for a vanishes there, hence so do all terms but the first one, 
and the first one vanishes as xd x q vanishes at x = 0. Thus, o~t, tm ([A, Q]) = ab + xe 
for some b G S£^( h T*X \ o), e G S , ^ m ( b T*A \ o). We deduce that there exists 
B G *™ _1 (A), E G * b n (A), i? G -V™ (X) such that [Q,A] = B(xD x ) + Ex + R. 
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As one can write E = E (xD x ) + J2 EjD Vj + R' with Ej,R' £ *™ _1 (A), and as 
x(xD x ),xD Vj £ Vo(A), the second claim is proved. □ 

Lemma 4.7. Suppose m > is an integer. Any A £ & b (X) defines a continuous 
linear map on H ' (X) by extension from C°°(X). 

Proof. We can use any collection B^ % ' £ Diff™(X), i — 1, . . . , N, such that at each 
point of °S*X at least one of the i?W is elliptic, to put a norm on H™' (X): 

WK^'ix) = E \\x- l B®u\\h W + WxMUx)- 

i 

We need to show then that for A as above, || Au\\ H m,i^ x ^ < G\\u\\ H m,i^ x y Since A 

is bounded on x~~ l L 2 (X), we only need to prove that for each i, \\x~ l B^' Au\\ < 
C'||w|| H m,i (x) . But x~ l B^A = Y J A j x- l B j with Aj £ and B 3 £ DifF^(X) 

by Lemma 1431 so \\x~ l B^Au\\ < Y,Cj\\ x ~ lB j u \\ as Aj are bounded on L 2 (X). 
This proves the corollary. □ 

As we work relative to x 1 Hq(X) = Hq (X), for k > we use the Sobolev spaces 

x l H^(X) = {u£ x l H£(X) : WA £ 1^(A), Au £ x l H^(X)}. 

These can be normed by taking any elliptic A £ 9^(X) and letting 

ll u IUff*; r (X) = W U L'H-(X) + \\ Au \\l'W (X)- 

Although the norm depends on the choice of A, different choices give equivalent 
norms. Indeed, if A £ ^(X), then let G £ *^ fc (A) be a parametrix for A, so 
GA = Id+E, AG = Id+F, E,F £ ^°°(X), and note that 

(4 x) \\M\x'h S (x) < \\AGAu\\ x i H r {x) + \\AEu\\ x i H r (x) 

< C(\\Au\\ x , HS(x) + \\u\\ x i HE{X )), 

where we used that AG £ ^°(X) and AE £ * b °°(X) C *b(A) are bounded on 
x 1 Hq(X) by Lemma l4~7l If A is elliptic, there is a similar estimate with the role of 
A and A interchanged, which shows the claimed equivalence. 

Lemma 4.8. If Q £ *"(A), then Q is bounded on x l H^(X). 

Proof. As Q is bounded on x 1 Hq(X), we only need to prove that for A £ ^(A), 
\\AQu\\ x i H5 (x) < C(\\u\\ x i H r {x) + \\Au\\ x , H r {x) ). But A = AQ £ <&£(X), though 
not necessarily elliptic, so by (|4.1[) . this estimate holds. □ 

Lemma 4.9. If L £ Difff(A) is elliptic, u £ x l H°'™(X), Lu £ x l H^(X), then 
u£x l H b s + k > oc (X). 

Proof. Let G £ ^^ k (X) be a parametrix for L so that GL = Id +R, R £ ^°°(X). 
Then u = G{Lu) - Ru. Now, if A € ^t{X) then Au = (AG){Lu) - (AR)u £ 
x l H s b ^{X) by Lemma El] since AG, AR £ fg(I), This proves the lemma. □ 

The conormal regularity theorem is global in each connected component of Y. 
It uses the following lemma, which shows that the boundary Laplacian commutes 
with P one order better (in terms of decay) than a priori expected: 
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Lemma 4.10. Let Ay € BiS 2 b (X) have normal operator given by Ay. Then 

[p, Ay] e xVifilmslix). 

Proof. Changing Ay by Q G £ Diff 2 , (AT) changes the commutator by an element of 
x Diffp Diff 2 (X) due to Lemma so the statement only depends on the normal 
operator of Ay. Similarly, it only depends on the normal operator of P. Thus, we 
may work on the model space [0, e) x x Y, replace P by (xD x ) 2 + i(n — \){xD x ) — 
x 2 A Y , Ay by Ay, and then the result is immediate. □ 

Proposition 4.11. Suppose I 6 R, u G x l H Q °°(X), Pu G C°°(X) and u G 
C°°(X°). Then for all e > 0, u G x l - e H^°{X) = x l - e H^°°(X). 

Remark 4.12. The proposition states that once one knows that u is smooth in X° 
and is in some weighted L 2 -space, one gets b-regularity relative to that space. 

Also, the proposition can be restated in terms of the standard b-spaces: u G 
x l +— £ H£°(X). The shift in the exponent is simply due to H£(X) being 
defined relative to Li(X), the P 2 -space relative to a non- vanishing b-measure. 

Proof. Assume first that I < -/(A). We prove that u G x l ~ e H^°°(X). We first 
note that by CorollaryEH u G H^°'(X) for all e > 0, i.e. we have full O-regularity. 
Let Ay be as above. 

As u G H™' l ~ e (X), Ay G x- 2 Diff^(A), we see that A Y u G H™^ 2 ^{X). Then 

(4.2) PAyw = K Y Pu + [P, A Y }u G F£°'' -1-£ (X) 

since [P, Ay] G a; Diff J Diff^X) C x^ 1 Diffg(X). (In fact, this can be phrased 
by saying that N(A Y ) and N(P) commute.) Thus, by Proposition 13.51 Ayu G 
H ™J-i-e( X y As [ xDx f u e H™' l - c (X), {{xD x ) 2 +A Y )u G H^'^-^X). Since 
(xD x ) 2 + Ay is elliptic in Diff^(X), Lemma B~l shows that u G x 1 ' 1 -^^ {X). 

Thus, 03 and [P,A Y ] G x Diff J Diff I (X) gives PAyu G x l ' e H^{X), so by 
Proposition 13.51 Ayw G H^°' ~ e (X). Proceeding as above, we deduce that u G 
x l -*Hl™{X). 

We now iterate this argument for A Y u in place of A Y u. So suppose we already 
know that u G x l ~ e H^ ~ lh °°(X) for all e > 0. Then [P,A Y ] G x Diff J Diff f C 

x- 1 DiflEgDiff2 ( *~ 1) (X), so 

PA Y u = A Y Pu + [P, A Y ]u G H^^^iX) 

Again, by Proposition SH A y u G H^^^^X). As (xD x ) 2k u G H™' l ~ £ {X), 
{{xD x ) 2k + A Y )u G ff£° ,i-:l-e pr). Using Lemma EH we conclude that u G 
aJ-i-efl-».~(jf). 

Equipped with this additional knowledge, we deduce that [P, A Y ]u G i?o°' ~ e (X), 
hence PAyit is in the same space. Applying Proposition 13.51 we see that Ayu G 
H^° ,l ~ e (X). Proceeding as above, we deduce that u G x l -*Hl^°°{X). This proves 
the proposition if I < —l(\). 

In general, if / > — £(A), we may apply the previous argument with I replaced 
by any V < -/(A) to conclude that u G x 1 ' H£°(X) for all I' < -/(A). Since 
u G x l L 2 (X), interpolation gives u G x L ~ e H£°(X) as stated. □ 
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We now consider P = □ — A acting on polyhomogcncous functions, or more 
generally symbols. Recall that u G A k (X) means that Lu G x L%(X) for all 
L G DiS h (X), so in particular u G x k L\[X). 

We remark that if s + , s_ 6 C with s + — s_ ^ Z, and a function w has the form 
x s +v + +x s -v^, v± G C°°(X), then the leading terms v±\Y (in fact, the full Taylor 
series of v±) is well-defined. However, if s + — s_ is an integer, this is no longer true, 
which explains some of the complications we face in stating the converse direction 
of the following lemma. 

Lemma 4.13. Suppose A G R, A ^ ( "~ 1)2 . Let 

be the (not necessarily real) indicial roots of (xD x +i(n— 1)){xD x )— A. //it G A k (X) 
for some k and Pu G C°°(X) and s+(A) — s_(A) is not an integer then there exists 
v± G C°°(X), sucft tftaf 

7/s+(A) — s_ (A) is an integer (in which case both s±(A) are rea(j t/ien t/ie analogous 
statement holds with v_ G C°°(X) replaced by 

w_ G C°°(X) +x s +( A )- s -( A )loga;C 00 (X). 

In either case, if v±\y vanish, then u G C°°(A'). 
Conversely, given g+,g~ G C°°(y) 7 tftere exist 

(i) u± G C°°(X) if s + (A) — s_(A) is not an integer, 

(ii) 

s+ (A)-s_(A)-l 

w+ G C°°(X), «_ - ^ a^' G x s +^- s -( A ) logxC 00 ^), 0j G C°°(Y), 

i/s + (A) — s_(A) is an integer, 
such that 

u = x s+w v + +x s - {x) v-, v±\ Y = g±, 

satisfies Pu G C°°(X). 

Proof. We start with the converse direction. As P = (xD x +i(n— 1)){xD x ) — X+Q, 
Q G xVmliX), for v G C°°{X), 

(4.3) P(x s v) = (s(n - 1 - s) - \)x s v + to, w G a;'' +1 C 00 (X). 

Thus, when s is an indicial root, P(x s v) G a; s+1 C°°(X) automatically, and otherwise 
given / G x s C°°(X), P(x s v) = f can be solved uniquely, modulo x s+1 C°°(X), with 
u G C°°(X). Iterating this argument, and using Borel summation, we deduce that 
unless the two indicial roots differ by an integer, given g+,g~ G C°°(Y), there exists 
v+,v- G C°°(X) such that 

u = x s+{x) v + +x s - (A) v_, v±\ Y =g±, 

satisfies G C°°(X). 
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If the two indicial roots differ by an integer (but are distinct, i.e. not equal to 
only a minor modification is needed in that we need to allow logarithmic 
factors. Thus, for v £ C°°(X), 

P(x s \ogxv) =(s(n — 1 — s) — A) \ogxx s v + (n — 1 — 2s)x s v + w, 

w £x s+1 \ogxC 00 (X) + x s+1 C 00 (X), 

so if s = s±(A), Pu = f,fe x s C°°[X), has a solution modulo x s+1 logxC°°(X) + 
x s+1 C oc (X), of the form u G x s \ogxC°°(X), so applying this with s = s+(A), the 
error term arising from s_(A) of the form a; s times a smooth function, can be solved 
away to leading order. Moreover, for s ^ s±(A), Pu — f, f S x s logxC°°(X) has a 
solution, modulo x s+1 \ogxC°°(X) + x s+1 C°°(X), of the form u £ x s \ogxC°°(X), 
so again iteration gives infinite order solvability in this case of the form: given 
eC°°(Y), there exists v + £ C°°(X), v_ £ C QC (X) + x s +^- s -^ \ogxC°°(X) 
such that 

u = x s+w v+ +x s - w v-, v±\ Y = g±, 

satisfies Pu £ C°°(X). 

On the other hand, suppose that u £ A k (X) and Pu £ C°°(X). As Qu £ 
A k+1 (X), we have ((xD x + i(n - l)){xD x ) - \)u £ A k+1 . Since near Y, using 
an product decomposition of a neighborhood of Y, A r (X) can be identified with 
C°°(Y;A r ([0,e))), we can treat Y as a parameter and solve this ODE. If there 
is no indicial root in (k, k + 1], one deduces that u £ A k+1 (X); otherwise u = 
Y]j x s j Qj + u' where the Sj are the indicial roots in the interval, gj are smooth and 
u' £ A k+1 . By the first part of the proof one can choose Vj as in the statement of 
the lemma (denoted by v± there) to get Uj = x Sj Vj £ A k with Puj £ C°°(X) and 
uj - x a ' g £ A k+1 . Thus, u-J2uj G A k+1 with P{u - J2 Uj) £ C°°(X), so one can 
proceed iteratively to finish the existence argument. Note that if gj\Y vanish, one 
concludes u £ A k+1 , which by iteration gives the uniqueness. □ 

In fact, the same argument also deals with the case A = (n — l) 2 /4, but as the 
result is of a slightly different form, we state it separately: 

Lemma 4.14. Suppose A = so s±(A) = -^rp-. If u £ A k {X) for some k 

and Pu £ C°°{X) then there exists v± £ C°°(X), such that 

u = x s+w v + +x s - {x) log XV-. 

Conversely, given g+,g~ £ C°°(Y), there exists v± £ C°°{X), such that 

u = x s +^v + + x 3 -^ \ogxv_, v±\y = g±, 

satisfies Pu £ C°°(X). 

Proof, s — s± (A) = (n — 1)/2 now satisfies s(n — 1 — s) — A — as n — 1 — 2s = 0, so 
(|4~3|) and imply that P(x s v 1 + x s logxv 2 ) £ x s+1 C°°{X) + x s+1 logxC°°(X). 
The argument of the previous lemma then shows the second claim. 

For the first claim, we need to observe that if u £ A k {X) and Pu £ C°°(X) then 
Qu £ A k+1 {X), so ((xD x +i(n-l))(xD x )-X)u £ A k+1 , i.e. {xD x +i(n~l)/2) 2 u £ 
A + . Proceeding as above, the only difference is that if s = £ (k, k + 1], one 
deduces that u = x s gi + x s log x gi + u', gj smooth, u' £ A k+1 . One finishes the 
proof exactly as above. □ 
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Since we already know (by virtue of Proposition 13.41 and Remark 13. 7|) that we 
can solve Pu' = f, f G C°°(X), with u' G C°°(X), modulo C^(X°), we deduce that 
these u can be further extended to be exact solutions near dX. 

5. Global solvability 

For global solvability, i.e. solvability on all of X rather than just near dX, of 
Pu = we need the additional assumptions (A1)-(A2). We thus assume that Y = 
Y + UY_, where Y± are unions of connected components of Y, and this decomposition 
satisfies that all bicharacteristics t i— > j(t) of P (i.e. those of □, independent of A) 
satisfy lim t ^ +oc j(t) G Y+, lim t ^_ 00 7(t) G Y_, or vice versa. In this case, noting 
that the sign of the \' term agrees with the others if r < min(0, 1 — 2Z(A)) (for they 
are all negative; recall /(A) = for the wave operator itself), one can easily 'cut 
and paste' the estimates with 

• near YL, r — r + > 1 + 2Z(A) (or just r = r + > max(0, 1 — 2Z(A)), r+ ^ 
1 + 21(A)), 

• near Y_, r = r_ < min(0, 1 — 2Z(A)) , and 

• standard microlocal propagation estimates in the interior of X 

to deduce that for a partition of unity X+ + X-+X0 = 1 with x+ supported near YL, 
identically 1 in a smaller neighborhood of YL, analogously with \-, Xo G C£°(A"°), 
there exists xo G C™(1°) such that 

(5.1) | |a; (r + -l)/2 x+1; ||2 + || a .(r-l)/2 x _ t;|| 3 + ^2 < C (||xo«|l^/ 2 + ll^f)- 

() () () o 

Let H™' q+ ' q ~ (X) be the space xV~ x q S H™(X), where x± are defining functions of 
Y±, we can put the norm 

on it. (Note that it is the completion of C co (X) with respect to this norm.) This 
is just x^H^(X) near Y±, H m {X°) in the interior. Let l± = (r± - l)/2. The 
argument of [51 Proof of Theorem 26.1.7] shows the following: 

Proposition 5.1. Suppose that A G BL 1+ > max(i, Z(A)), /_ < — max(i, Z(A)). 
Then 

N l+t i_ ={ve H^'- l +'- l -{X) : Pv = 0} 

is finite dimensional, and for f G H ' + ' _ (X), f orthogonal to Ni +> i_ , Pu = f has 
a solution u G Hq' 1+ ' 1 ~ (X). 

Moreover, elements of Ni + i_ are in H™' l ~ l ~(X) for all I < -l+, are Schwartz 
at YL , and have an expansion as in Lemma \4-13\ at Y + . 

Remark 5.2. Note that the expansion of Lemma [4.131 implies that Ni + j_ are in 
H^' l '°°(X) for all I < -Z(A), not merely I < 

Proof. We first prove the last statement. For v G Ni + j_, by Corollary 13.61 v is 
Schwartz at Y_. In particular, v is C°° near Y_, so by the standard propagation 
of singularities for P, v G C°°(A°). Then, by Corollary EH v G H™' l '~ l ~ for all 
I < — Z+. By Proposition 14.111 and the remark following it, u G x l+ ^~ H£°(X) = 
A l+S ^-(X) for all / < -l+. Thus, by Lemma ECU it has an expansion at Y+ of 
the form given by Lemma fl.lSI 
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This in particular implies that the commutator calculations giving rise to (|5.1j) 
can be applied directly (without mollification) to all v E Ni +t i_ The proof of the 
first part is finished as in [5] , and the second part can then be proved exactly as in 

i- □ 

Note that the role of Y± is reversible, so the estimates, hence the proposition, 
also hold with l± interchanged. Correspondingly, we deduce that the solution u of 
Pu = f above is unique modulo the finite dimensional space A r _; + _;_ . 

One can also get uniqueness, namely that 

Proposition 5.3. Suppose u G C°°(X) and Pu = 0. Then u = 0. 
In fact, it suffices to assume that u is Schwartz at Y + . 

If we merely assume that u is Schwartz at a connected component Yj of Y , and 
Pu = near Yj , then we can still conclude that u = near Yj . 

Proof. The proof is very similar to |T?J Section 4] and to [TH]- Consider Ph = 
x~ 1 ' h h 2 Px 1 ' h . The basic claim is that the semiclassical symbols of ReP/, 6 
Diffp h (X) and ImP/> 6 Diff^ (X) never vanish at the same place at Y . In fact, as 
P is formally self-adjoint, one has 

P h = h 2 P + x- 1 ' h [h 2 P.x 1 /% 
ReP h = h 2 P +^[x- 1 /\[h 2 P,x 1 / h }}, 

lmP h = ^{x- 1/h [h 2 P,x^ h } + [h 2 P,x 1/h ]x- 1 l h ). 

Now, for Q 6 x l Diff^(X), x-^ h [Q,x^ h ] e x l Diff^X), so if we only want to 
compute the commutators modulo higher order terms in x, we can work with the 
normal operator of P instead of P. Also, modulo higher order terms in h, only the 
principal symbol of P matters in the calculations, as we are considering h 2 P, and 
changing P by a first order term changes h 2 P by an element of h Diff q ^(X). Thus, 
a straightforward computation gives 

ReP^ = (hxD x ) 2 - x 2 A Y + \[x- 1,h , [h 2 (xD x ) 2 ,x 1 / h ]} + R, 

= (hxD x ) 2 - h 2 x 2 A Y - 1 + Pi, 
ImP, = ^(x-y h [(hxD x ) 2 ,x^ h ] + [{hxD x ) 2 ,x 1 /' l }x- 1 / h ) + R 2 = -2hxD x + R 2 , 

with Pi G hI)iS 2 h {X) + xDiS 2 h (X), R 2 G hDiSl h (X) + xBiS^ h {X). Moreover, 

i[ReP h .lmP h ] = i[-h 2 x 2 A Y , -2hxD x ] + P 3 = -4h 3 x 2 A Y + hR 3 , 

P 3 £ hDiS 2 h (X) + xDitt 2 h (X). Thus, 

z[Re P h , Im =h + 4hReP h - h(lm P h ) 2 + hR 4 , 

with P4 having the same properties as P3. 

Now let u h = x- x l h u G C°°(X), so P h u h = and 

0= \\P h u h \\ 2 = HReP^fcl^ + IIImP^f + ^tRePh^PhK,^) 

= || ReP h u h f + (1 - h)\\ lmP h u h \\ 2 + h\\u h \\ 2 + 4h{Re P h u h , u h ) + h(R 4 u h ,u h ). 

This is the analogue of Equations (4.2) and (4.3) of [17], except that here terms 
arising from the commutator i[Re Ph,Im Ph] do not have an additional factor of x 
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compared to the first two squares on the right hand side. The proof can be finished 
exactly as in [17], writing R4 = hRs + x 1 / 2 R^x 1 / 2 , R5, Rq € Di& h (X), and noting 
that — RePfj + (Imf\) 2 is elliptic second order, so 

\(hR 5 u h ,u h )\ < Ch(\\ ReP h u h \\ \\u h \\ + \\ lvaP h u h \\ 2 + \\u h f), 

\(x^ 2 R 6 x^ 2 u h ,u h )\ 

< Ch{\\ RePhX^UhW \\x^ 2 u h \\ + || lmP h x l ' 2 u h \\ 2 + \\x^ 2 u h \\ 2 ) 

< C'h(\\ ReP h u h \\ Wx^uhW + \\lmP h u h \\ 2 + \\x^ 2 u h \\ 2 ). 
Indeed, for S > one writes 

\\xV>u h f = \\ X 1/2 u h \\l<s + \\x 1/2 u h \\l> s < S\\u h \\ 2 + S'- 2 / h \\uf, 

so 

> (1 - dh)\\ ReP h u\\ 2 + (1 - C 2 h)\\ lmP h u\\ 2 + h{l - C 3 h - C 4 S)\\u h \\ 2 

-C^- 2 / h \\u\\ 2 . 

Thus, there exists ho > such that for h G (0, ho), 

hC 5 5 1 - 2 l h \\u\\ 2 >h{ l --C i 5)\\u h \\ 2 . 

Suppose 8 € (0,min(^-, -£-)) and suppwnja; < |} is non-empty. Then ||w/i|| 2 > 
C 6 (5/4)~ 2 / h with C* 6 > O.'Thus, 

C 5 S\\u\\ 2 > ^4 2 /'\ 

As the right hand side goes to +00 as h — ► 0, this provides a contradiction. 

Thus, u vanishes for x < 8 /A, and then the usual hyperbolic uniqueness (well- 
posedness of the non-characteristic Cauchy problem) gives that it vanishes on X. 

□ 

Combined with Proposition 15. II this gives: 

Theorem 5.4. Suppose that A G M, l + > max(|,Z(A)), Z_ < — max(i, Z(A)). Then 
for f G Hq' 1+ ' 1 ~ (X), Pu — f has a unique solution u G Hq' 1+ ' 1 ~ (X). 

Proof. With the notation of Proposition ^. li we want to prove Ni+ i_ = {0}. But 
for v G Ni +t i_, by Corollary 13. 6\ v is Schwartz at Y-. Thus, by Proposition 15. 3\ 
v = 0. Thus, by Proposition l5.il the required u exists. 

Conversely, if u G Hq' 1+ ' 1 ~ (X) and Pu = then by Corollary 13. 61 u is Schwartz 
at Y + , so by Proposition [531 u = 0. □ 

We also deduce: 

Theorem 5.5. Suppose A ^ . Given g± G C°°{Y + ) there exists a unique 

u G C°°(X°) such that Pu = and which is of the form 

u = x s +^v+ + x s -( A V, v ± \ Y+ = g ± , v+ G C°°(X), 

s+ (A)-s_(A)-l 

v-- a 3 x3 e ^ +(A)_S - (A) logxC°°(X), a 3 G C°°(Y±). 

j=o 

If S-|-(A) — S-(A) is not an integer, then w_ G C°°(X). 
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r j \ 2 

On the other hand, if X = ~ , then given g± G C°°(Y + ) there exists a unique 
u G C°°(X°) such that Pu = and which is of the form 

u = x in - 1)/2 v + + x [n - 1)/2 logxv_, v±\ Y+ =g±, v±eC 00 (X). 

/ ^\2 

Proof. Suppose A / (n ~. ' . As shown in Lemma [4.131 there exists uq supported 
near Y + and of the desired form there, such that Puq € C°°(X). By Theorem 15. 4[ 
for any l + > max(i,Z(A)) and Z_ < — max(i,^(A)) there exists a unique u\ G 
Hq' 1+ ' 1 {X) such that Pui = —Pu G C° C (X). As l± are arbitrary subject to the 
constraints, and u\ is unique, ui G (X) /or 1+ > max(±,Z(A)) Z_ < 

-max(|,i(A)). By Corollary EH Mi is Schwartz at Y" + . Thus, u = uq + u\ satisfies 
Pu = 0, and is smooth near Y+, so by the standard propagation of singularities 
u G C°°{X°). As u G ffo'^PO for all Z_ < - max(±, Z(A)) near YL, CorollaryO 
gives u G H^°' l ~ (X) for all such By Proposition 14. 1 II and the remark following 
it, u G a; i - + T i iJ fc 00 (A) = A l -+^(X) for all such L_. Thus, by Lemma [TT3] it 
has the stated form near F_ . 

Conversely, if u has the stated properties and g± = 0, then v± are Schwartz at 
Y + by Lemma T4. 131 so u is Schwartz at Y+. Then u = by Proposition [531 

If A = the same argument, but using Lemma [4 . 1 41 instead of Lemma l4.131 

completes the proof of the theorem. □ 

6. The Cauchy problem 

We now consider global solutions for the Cauchy problem posed near Y±. 

Let T be a compactified time function, as in the introduction. For any constant 
to G (—1, 1), and a vector field V transversal to St 0! P is strictly hyperbolic, and 
the Cauchy problem 

Pu = in X°, 

(6.1) u| Sto =^ , 

Vu\ StQ = 

^0,^1 € C°°(S to ) is well posed. 

Theorem 6.1. Let s±(A) = 2=i ± y / ( "~ 1)2 - A. Assuming (Al) and (A2), the 
solution u of the Cauchy problem (|6.ip has the form 

(6.2) u = x s + w v + +x s - {x) v-, v ± eC°°{X), 

if s_|_(A) — s-(A) = 2^/ ( -"~ 1 - ) A is not an integer. If s+(A) — s_(A) is an in- 
teger, t/ie same conclusion holds if we replace V- G C°°(X) by v- G C°°(X) + 

2;S+ (A)- S _(A) \ ogx C°°(X). 

Proof. As P is strictly hyperbolic with respect to S to , [HI Theorem 23.2.4] guar- 
antees the existence of u G C^{X°) with Pu = in a neighborhood of 5t and 
having the required Cauchy data. We may choose ti < t < t 2 so that Pu a = 
for T G (ti,i2)- Let Xi>X2 £ C°°(X) be such that xi = 1 m a neighborhood of 
T > to, xi is supported in T > ti, while X2 = 1 hi a neighborhood of T < to, 
supported in T < t 2 . In particular, X1X2 is supported where T G (ti,ia), is iden- 
tically 1 near St , and each Xi is identically 1 on the support of the dXj, j 7^ i- 
Then P(xiX2«o) = [P, Xi] w o + [P X2]"o- Denoting these two terms by /1, resp. f 2 , 
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we use Theorem 15.41 to solve away f\ towards Y+ and /2 towards Y- so that the 
Cauchy data are unchanged. 

First, by Theorem 15.41 with any I > max(|, /(A)), there exists U2 G Hq' '~ (X) 
such that Pu2 = fs- By Corollary 13.61 u^ is Schwartz at Y + , and then by Propo- 
sition [33J u 2 = near Y + . Hyperbolic propagation then shows that supp it 2 C 
{T > to} as f% is supported in this set, so ui = near St . In addition, as in 
the argument of Theorem [53] we deduce that u-i G C°°(X°) has an expansion as in 
Theorem 15.51 

Interchanging the weights at Y±, we can similarly show the existence of u% G 
H ' ' (X) such that Pu\ = /i, suppui C {T < to}, and m having an expansion 
at Y+. Thus, u = X1X2U0 — u\ — u-i € C°°(X°) satisfies Pu — 0, u|s t = 4>o, 
Vu\s t = ipii and u has an asymptotic expansion as in Theorem 15.51 proving the 
existence part. 

Uniqueness follows easily, for if u solves the Cauchy problem with ^0 = 0, = 0, 
then u = near S to , hence vanishes globally. □ 

It is useful to relate the Cauchy data at different hypersurfaces to each other, 
particularly for hypersurfaces near Y +1 resp, Y-. This is very easy using the stan- 
dard FIO result. We renormalize this operator in order to make all entries in the 
FIO matrix have the same order. Namely, let A t . be the Laplacian of the restriction 
of g to S tj , j = 1,2, so A tj > as S tj is space like. Let A' t . denote the operator 
which is At j on the orthocomplement of the nullspace of A^. and is the identity on 
the nullspace, so A' tj is positive and invertible. 

Proposition 6.2. For any t\,t% G (—1,1), the map C' tl ,t 2 sending Cauchy 

data of global smooth solutions of Pu = at St t to Cauchy data at St 2 : 

C tut2 : {{A' ti )^u\ SH ,Vu\ Stx ) ~ ((A^uk^ukJ 

is an invertible Fourier integral operator of order corresponding to the bicharac- 
teristic flow. 

7. The scattering operator 

In order to prove that the scattering operator is a Fourier integral operator, we 
construct a parametrix as a conormal distribution on a resolution of X x Y + for the 
solution operator, also called the 'Poisson operator', gJ) i— > u with notation as 
in (JTTTJ) and (TO]) , 

Near Y+, this can be done by considering [X x Y" + ;diag y+ ]. On this space the 
parametrix is a conormal distribution near Y+ associated to the 'fiowout' of points 
in Y + . That is, for q' G Y + , consider the bicharacteristics approaching °S*,X. 
These form a Lagrangian submanifold of T*X°, which near Y + has constant rank 
projection (since the rank at the front face is maximal, namely n — 1), and is thus 
the conormal bundle of a submanifold F g i of X. These F q i depend smoothly on 
q' so that F = Uq'Fq' x {q'} is a smooth submanifold of X° x Y+, and indeed it 
extends to be smooth to [X x F + ;diag y+ ]. 

In order to orient ourselves, we first make some remarks regarding distributions 
conormal to F. First, recall that if M is a manifold with corners of dimension 
to, and Z is an interior p-submanifold, I P (M,Z) is the space of distributions on 
M conormal to Z, see [TU1 E]. Here we only need the case where Z meets all 
boundary faces transversally; in fact, in this case, Z only meets a (codimension 
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one) boundary hypersurface. Thus, in local coordinates (x,y), x — (xi, . . . , Xk), 
y = (j/i, . . . , y m -k) in which M is locally given by Xj > for all j, and Z is given 
by Vi = ■ ■ ■ = Vn — 0, elements of I P (M, Z) have the form 



with a e 5 p+ ^ 2JV )/ 4 (M;R JV ), j/' = (yi, . . . , y N ). Note that x behaves as a param- 
eter, i.e. the presence of boundaries does not cause any complications, hence the 
standard treatment in the boundaryless case [7J [5] actually suffices. Note that if 
A G DifT(Af) and u G I P (M, Z) then Au G I p+m (M, Z), and if 4 is characteristic 
in Z, i.e. its principal symbol vanishes on N*Z, then Au G I p+m ~ 1 (M, Z), with 
CTp-|- m — i (j4u) = H a a m (u) + 6u, where 6 depends on A only. This equation is an 
ODE along the bicharacteristics of A, and is called a transport equation. 

We also need to allow weights, i.e. consider the spaces x s I p (M,Z). Diff (M;Z) 
is not well-behaved on these spaces (because of derivatives possibly falling on x s ) 
but Diff b (A/) is. 

Lemma 7.1. (see ^[Section 18.2] and [10 ) Suppose that A G Diff™(M). Then 

(7.1) A:x s I p (M,Z)^x s I p+m (M,Z). 
If A is characteristic on Z , then 

(7.2) A : x s I p (M, Z) — au'J^-^Af, Z), 

and there is function b depending on A only such that a p+m ^i(Au) = H a a m (u)+bu. 

Proof. As x~ s Ax s G Diff b "(Af) C Diff m (Af), CLX) follows immediately from the 
remarks above. Next, if A is characteristic on Z, then so is x^Ae", so the remarks 
above prove (|T. 2[) . As the principal symbol of x~ s Ax s is the same as that of A, 
a p + m -i(Au) = H a o~ m (u) + bu follows. □ 

In our case, M = [X xY + ; diag y+ ] , and Z = F. The transport equation will allow 
us to solve away errors modulo smooth terms in our construction of the 'Poisson 
operator', (g +1 gJ) i— > u. However, we need to see first what the 'errors' are errors 
of, i.e. where the Schwartz kernel of the Poisson operator comes from, which will 
also give a relationship between the orders s and p above. 

Even for arbitrary Y, the model on the front face is the same as when Y is 
Euclidean space with a translation-invariant metric. Let y denote local coordinates 
on Y, as well as their extension to X, so (x,y) are local coordinates on X. On 
X x Y + then we have local coordinates (x, y, y'), where y 1 is the pull-back of y from 
the second factor. (The pull-back of y from the first factor, X, is still denoted by 
y.) Using projective coordinates 

X = x, Y = , y , 

x 

the P = □ — A becomes 

(XD X - YD Y +i(n- l))(XD x - YD Y ) - ^ h t ,(y')D Y% D Y] - A, 

modulo XDiff b ([X x Y+;diag y+ ]). To analyze this operator for fixed y' , we may 
arrange that hij(y') = Sij, so the operator becomes 

(7.3) (XD X - YD Y + i{n - 1))(XD X - YD Y ) - Ay - A. 
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When acting on functions of the form u = x s v, v a function of Y, XDx becomes 
a multiplication operator, and the operator we arrive at after this substitution is 
a degenerate PDE with radial points over \Y\ = 1, i.e. where F hits the front 
face. This is indeed what enables us to find solutions supported in \Y\ < 1, with 
singularities carried away by F. 

While this form is helpful in seeing the big picture, we need to solve this exactly 
at X — to leading order, for which it is useful to view □ on the warped product 
model as the analytic continuation of the Laplacian on hyperbolic space, which 
is arrived at by complex rotation in x (replacing x by ix), i.e. considering the 
Laplacian of dx J[ h . Correspondingly, the explicit solutions we are interested in are 
analytic continuations of the Eisenstein functions (Poisson kernel) on hyperbolic 
space, i.e. they take the form 

X S (\Y\ 2 - 1 ± i0) s , -s(n + s - 1) = A. 

Note that these values of s arc different from the usual indicial roots; these give 

a = s±(A) = ± \Z(^) 2 - A = S ±( A ) - (» - !)■ 

We in fact have two interesting solutions corresponding to branches of the analytic 
continuation. As we are interested in solutions supported inside \Y\ < 1, we take 
their difference, 

x s [{\y\ 2 - 1 + ioy - (\y\ 2 - 1 - ioy] = Cs x s (\y\ 2 - iy_, 

with c s = e l7TS — e~ l7TS if s is not a negative integer, and 

X S [(\Y\ 2 - 1 + iOy - (\Y\ 2 - 1 - ioy] = cX'S^'-^QYl 2 - 1), 

with c s = ^rr^rjji if s is a negative integer. Here the notation is that if / is a 
distribution on R which is conormal to the origin, then /(|F| 2 — 1) denotes T*f, 
where T : M™ _1 — > R is the map T(Y) = \Y\ 2 — 1. The prcimagc of the origin under 
T is the unit sphere, and on the unit sphere the differential of T is surjective, so 
the pull-back of these conormal distributions indeed makes sense. 

If the boundary is actually Euclidean, then near Y + x Y + we thus obtain an exact 
solution with singularities on F, 

£o,±(z,2/,y', A) - C s x s (l - ljL ^f-) s + , 
with C s to be determined and s = s±(A), if s is not a negative integer, and 

Eo,±(x,y, y', A) = CWS^ (1 - ^|^) 
if s is a negative integer. Note that for each A, 

E ,± = £ ,±(A) e x s I m ^([X x y+;diag y+ ],F), m(s) = -s - ^±1 s - s±(A). 

Lemma 7.2. Suppose that s±(A) ^ — — N+. Then there is a constant C s ^ 
such that for all <fi € C°°(Y + ) the operator E ±(A) with Schwartz kernel E± alh: 

Eo,±^= f E , ± (x,y,y',\)cb(y')dh(y') 
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satisfies 

Eo,±4> = x s[X) v, v G C°°(X), v\ Y+ = <t>. 

Remark 7.3. Note that for A > ^£ - 1, the condition s±(A) ^ - N+ 

automatically holds. For □ itself (i.e. A = 0) the condition holds if n is even. In 
addition, the condition always holds for one of the two indicial roots, namely the 
larger one (i.e. the one with more decay/less growth at Y+). 

Proof. Suppose first that s±(A) is not a negative integer. 

Changing variables in the integral we deduce that for <fi e C£°(Y+), and s = s±(A) 
still, 

J E 0!± (x,y,y',\)<t>(y')dy' = x n ~ 1+8 ^ J (1 - \Y\ 2 ) s + ^y - xY) dY 

= x s ^v, v e C°°(X), v(0, y) = C s ((l - \Y\ 2 Y +1 l)0(y), 

where the second factor in the expression for v(0, y) is the evaluation of the distri- 
bution (1 — |F| 2 )^_ on 1, and where we used that s±(A) = s±(A) + (n — 1). We need 
to check for which values of s does C s vanish, so we compute this pairing. 

For Res > —1, the distributional pairing is an absolutely convergent integral, 
which in polar coordinates becomes 

c„_2„,n-i , „ c„_ 2 r(^)r( s + i) 



Cn-2 J (I- P 2 )^"- 1 dp = C ^1 B (^1,S + 1) 



2r(^ + s + i) 



where c„_2 is the volume of the (n — 2)-sphere and B is the beta-function. As both 
the distributional pairing and the T function are meromorphic in s (indeed analytic 
away from — N), we deduce that 

((1 l y U+'^- 2 r(V + * + i) 

for all s which are not negative integers. This vanishes only if s e — — N + and 
n is even (so s is not a negative integer). 

If s = s±(A) is a negative integer, say s = — k, 

J E , ± (x,y,y',X)cP(y')dy' =x n ~ 1+s ±W J (1 - \Y\ 2 ) s + 4>(y - xY) dY 

= x s ^v, veC^(X), v(0,y) = C s (6 { Q - s - 1 \l-\Y\ 2 ),l)cf>( y ). 
The distributional pairing now becomes 

^(<T 1} (.), (1 - ,) ( - 3)/2 ) = ^(^T(l - -) ( - 3)/2 Uo. 
If n is even, all derivatives of (1 — z)^" -3 ^ 2 at z = are non-zero, while if n is 

2 



odd, the derivatives of order < 2— - are non-zero, so this pairing vanishes only if 



8 = -k € -N+. 

Combining these two cases, s = s±(A) ^ — — N + implies that the respective 
distributional pairings are non-zero. Letting C s to be their reciprocal yield E .±(X) 
satisfying the lemma. □ 

If the metric is not exact warped product, then £?o,± will play the role of the 
model at the front face of [X x Y+; diagy + ], which then will need to be 'extended' 
into the interior. First, let y : Y + x Y + — > be local coordinates on the first 
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factor of y+ centered at the diagonal so that at the diagonal, the metric h lifted 
from the first factor is the standard Euclidean metric dy 2 . That is, informally, 
y = y(y') is a family of local coordinates on Y + , parameterized by y' G !+, so 
that for fixed y', y(y') gives local coordinates centered at y' in which h is dy 2 at 
the center, y(y') = 0. Thus, with the notation considered above in the Euclidean 
setting, we can take y — y — y' . Let Y = — - , so (x, Y, y') form a local coordinate 
system in a neighborhood of the interior of the front face of [X x Y + ; diagy ]. 

As F is a C°° codimension 1 submanifold of [X x Y + ; diag y+ ] transversal to the 
front face, intersecting it in the sphere \Y\ = 1, there exists a C°° function p on 
[A" x F + ;diag y+ ] such that p defines F (i.e. p vanishes exactly on F, and dp does 
not vanish there), and p|g = 1 — \Y\ 2 . We let r > be defined by r = (1 — p) 1 ^ 2 , so 
r = \Y\ at ff, and for convenience we often write (slightly imprecisely) (1 — r 2 ) s + , 
etc., for p s + . Our model is then 

E 0t± (x,y,y',X) = C s x s (l-r 2 )% = C s x s p s + , 

if s = s±(A) is not a negative integer, and 

Eo,±(x,y,y', A) = C.z*^— x) (l - r 2 ) = C^S^^ip) 

if s is a negative integer, with C s as in Lemma 17.21 
Then we want to find 

E± g a; s 7 m W([X x r+;diag y+ ],F), m(s) = - 8 -*±±±, s = J±(A), 

with PE± =Q, E±- E 0y± G s s+1 7 m W([Jf x y + ; diagy ], F), and f? ± vanishing 
to infinite order off the front face. The equation PE± G C°°(A" x y + ) becomes a 
degenerate transport equation at the level of principal symbols and can be solved 
to leading order. In fact, in order to simplify the transport equation, which is an 
equation for the principal symbol of E±, given by an ODE along the Lagrangian, 
N*F, it is convenient to notice that we want 

E± = ax s (l-r 2 ) s + +E' ± , a £ C°°([X x r + ;diag y+ ]), 

E' ± G x s I m ^- 1+e {[X x y +; diag y+ ],F), 

e G (0, 1) arbitrarily small, so the principal symbol of E± can be identified with 
o,\f , and the transport equation is an ODE for a\p. Namely, 

PE ± = {Qa)x s {l -r 2 )*- 1 + E±, E± e x s I m ^+'([X x Y + ; diagy + ], F), 

where Q is a first order differential operator of the form Q = xV + b, V a vector 
field tangent to F transversal to OF - xV(q) is a non- vanishing multiple of the 
push-forward of the Hamilton vector field H p evaluated at the one-dimensional 
space N*F q i \0. (This vector field is homogeneous, so the choice of a G N*F q i only 
changes the push forward by a non- vanishing factor.) 

Solving the transport equation and iterating the construction gives a new E± G 
x s I m ( s )(JX x Y + ; diagy + ] , F) vanishing to infinite order off the front face with 
PE± G x s+1 C°°([X x F + ; diagy + ]); we show this in Proposition 17.81 below. In 
fact, we can do better: we can ensure that near Y + (where this makes sense) E± 
is supported in the interior of the light cone] this is important as we show momen- 
tarily. 
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In order to remove the leading term at the front face (i.e. to improve the error, 
PE±, to x s+2 C°°([X x K)_;diagy • ]), which can then be further iterated away), we 
need to study P acting on functions of the form x a v, v £ C°°([X x F + ; diagy- + ]), 
modulo x a+1 C°°{[X x Y + ; diagy- ]). This only uses the model at ff. But ([7T5]> gives 

x- a Px a v = P a v, P a = (YD Y -i(n-l- <r))(YD Y + ia) - Ay - A, 

with P(j on operator on Euclidean space identified with the fiber of the front face 
over y' . This is of course a differential operator with smooth coefficients, but it is not 
elliptic. To see its precise behavior, it is convenient to introduce polar coordinates 
(r, Lj) in Y. (This agrees with our preceeding definition of r at the front face.) In 
such coordinates, 

n — 2 1 

Pa = (rD r -i(n-l- a)){rD r + ia) - D 2 r + i — — D r — — A w — A, 

with Au the positive Laplacian on the standard (n — 2)-sphere. The principal 
symbol of P a is (r 2 — 1)|£| 2 — r ~ 2 M5, with (£,77) denoting the dual variables of 
(r, lu). Thus, P a is elliptic for r < 1, i.e. inside the light cone. A straightforward 
calculation shows that P a is microhyperbolic for r > 1; it has some radial points at 
r = 1. There are two slightly different (but related) aspects of P a to address: the 
solvability of the transport equations, i.e. the removability of singularities at r = 1, 
and the solvability of smooth terms. 

We start with the transport equations. It is convenient to consider the conjugate 
(1 — r 2 )~ s P CT (l — r 2 ) s , more precisely, in view of the singularity of the conjugating 
factor, (1 — r 2 ± iO) _s P (T (l — r 2 ± i0) s , considered on all of the front face, i.e. 
as an operator from C°°(ff) to C~°°(ff). The following lemma is the result of a 
straightforward calculation when replacing ±i0 by ±ie, and the lemma then follows 
by taking the limit. 

Lemma 7.4. For all s E K., P a satisfies 
(1 - r 2 ± iOy s P a (l - r 2 ± i0) s 

n — 1 

= 4s(s - cr)(l - r 2 ± iO)- 1 + ( P(j - 4s{rd r + s - a + — — )) 

= 4s(s -<j)(l-r 2 ± tO)- 1 + P a _ 2s 

as operators from C°°(ff) to C"°°(ff). 

We in fact always need logarithmic terms to solve away singularities because 
there are automatic integer coincidences between the powers of x we need in the 
Taylor series, i.e. a, and the orders of the singularities along F, i.e. s. 

Lemma 7.5. For all s el, k S N, P a satisfies 

Pail - r 2 ± i0) s log(l - r 2 ± i0) k 

=4fc(2s - cr)(l - r 2 ± iOy- 1 log(l - r 2 ± i0) fc " 1 

+ 4s(s - cr)(l - r 2 ± i0) s_1 log(l - r 2 ± i0) k 

(7.4) + (1 - r 2 ± i0) s log(l ~r 2 ± iO) fe P CT _ 2s 

fe-2 

+ (1 - r 2 ± i0) s log(l - r 2 ± zO)*-^^!, 
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as operators from C°°(ff) to C °°(ff), where the Qj, j = 0, . . . , k — 1 are first order 
differential operators with smooth coefficients on ff (depending smoothly on s, a, k). 

Remark 7.6. The principal utility of allowing logarithmic singularities arises if s — 
er, in which case the second term on the right hand side is missing, hence the first 
term can be used to remove error terms with a lower power of logarithm (that could 
not be removed without logarithms, i.e. by the preceeding lemma). 

Proof. The case k — follows from the preceeding lemma. We then proceed by 
induction. If k > 1, and the result has been proved for k replaced by k — 1, then 
for a e C°°(ff), 

PJl - r 2 ± i0) s log(l - r 2 ± i0) k a = -—PJl - r 2 ± i0) s logfl - r 2 ± i0)* _1 a 

ds 

shows that we simply need to differentiate (|7.4[) (with k — 1 in place of k) with 
respect to s. The only terms giving rise to additional factors of logarithms are the 
ones in which (1 — r 2 ± i0) s_1 or (1 — r 2 ± i0) s is differentiated. As we are applying 
the result with k replaced by k — 1, the last two (residual) terms of (|T.4[) (for k — 1) 
give rise to residual terms (for k). Also, the only term that is not negligible even 
though it has a power of logarithm less than k is the first one, with a factor of 
(1 — r 2 ± iOy^ 1 . Thus, the first three terms of (|7.4p (for k — 1) will contribute to 
the last two residual terms (for k) except when (1 — r 2 ± «0) s_1 or (1 — r 2 ± i0) s 
is differentiated, or when the coefficient of the second term is differentiated. The 
latter gives 4(2s — <r)(l — r 2 ± z0) s_1 log(l — r 2 ± zO)*^ 1 , so altogether we have 
4(2s - a) + 4(2s - a)(k - 1) = 4fc(2s - a) of (1 - r 2 ± iO) 3 " 1 log(l - r 2 ± i0)* _1 , 
giving the desired result. □ 

Corollary 7.7. If s ^ — 1 and s + 1 ^ a then for b smooth function on S™ -2 , there 
exists a unique smooth function q on §"~ 2 such that 

P(x a (l - r 2 ); +1 g ) = x°(l - r 2 ) s + (b + (1 - r 2 )e), 

holds near E> n ~ 2 , with e smooth near §™~ 2 . 

More generally, under the same assumptions on s, for b smooth function on 
E> n ~ 2 , there exists a unique smooth function q on S™ -2 such that 

P(x' T (l-r 2 );+ 1 log(l-r 2 )^) 

fc-i 

= x°(l - r 2 );(log(l - r 2 )\b + log(l - r 2 ) J + e 3 + (1 - r 2 ) log(l - r 2 )\e), 

holds near S™~ 2 , with e, ej smooth near S™ -2 , j = 0, 1, . . . , k — 1. 

If s = — 1 or s + 1 = a, but a y 0, £/ien /or 6 smooth function on S n ~ 2 , there 
exists a unique smooth function q on §"~ 2 such that 

P(^(l-r 2 ); +1 log(l-r 2 )^ 1 g ) 

k 

= X°(1- r 2 );(log(l - r 2 )\b + log(l - r 2 )\e 3 + (1 - r 2 ) log(l - r 2 )^ +1 e), 
ZioWs near §' l_2 ; with e, smooth near § n_2 ; j = 0, 1, . . . , k. 

Proof. In the first case, let q = 4s _1 (s — cr) _1 6, and apply Lemma T7.41 expressing 
(1 — r 2 )^ +1 as a difference of (1 — r 2 ± i0) s+1 . Uniqueness is clear. 

In the second case, proceed the same way, applying Lemma 1731 □ 
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Proposition 7.8. The transport equations can be solved near the front face, i.e. 
there exist 

E ± G x s I m ^([X x r+;diag y+ ],F), m(s) = -« - ?2±± s = s±(A), 

with 

(7.5) E ± -E . ± G z s+1 / m «([X x ^+;diag y+ ],F), 

P.E± G a; s+1 C°°([X x y + ; diagy ] ) , and J?± vanishing to infinite order off the front 
face. 

Proof. First, with any E± G x s / m ( s )(LY x F + ; diag y+ ] , F) extending P ,± in the 

sense of (|7.5|) . having an expansion in terms of (1— r 2 )+, so E = x s± ^ (1— r 2 )^ a, 
a smooth, a| x =o = 1, one has 

P£±=^( A ) +1 (l-r 2 ); ±(AH 5, 

b smooth. By the corollary (if s±(A) ^ 0), one can find E\± — x s± ^ +1 (l — 
r 2 )+ ±(A) 6 such that 

PE lt ± - PE± G ^±( A ' +1 (1 - r 2 )f W C°° + ^±( A )+ 2 (1 - r 2 )^^^, 

so replacing P± by E± — £a,±, one has an extension of Eo } ± of the same form as 
the original E± , but with 

PE± G ^±( A ) +1 (1 - r 2 ); ±(A) C°° + ^±( A )+ 2 (1 - r 2 )^ - ^ 00 . 

Leaving the first term unchanged, one iterates the second term away, using Ej t ± — 
x*±( A )+J(l-r 2 )+ ±(A) 6 to remove errors in x s ^ +j (l-r 2 )^ 1 ^' 1 ^ , with the result 
that the new E± satisfies 

PE ± G ^±( A ) +1 (1 - r% ±W C°° + x S ±M +j+1 {l - r 2 )^- 1 ^. 

Note that there is no obstacle for this procedure as long as s±(A) ^ 0. By an 
asymptotic summation argument one gets an E with 

PE ± G x s ^ +1 {l - r 2 ) S + ±W C°° + (1 - r 2 ) ^ ±{A) - 1 d~. 

For the last term the singular transport equations are now easily solvable, so one 
obtains near the front face 

PE±ex s ±M +1 (l-r 2 ) $ ± W C°°. 

Now using the corollary we can find E\ — x s± ^ A ' +1 (l — r 2^s±W+ 1 \ g(\ _ r 2 ) + 
such that 

PE h ± - PE ± e^±( A ) +1 (i - r 2 ) ;±( A )+ 1 C oo 

+ ^±( A )+i(l - r 2 )^ (A)+1 log(l - r 2 ) + C°° 

+ x ^W+ 2 (l-r 2 Y + ±W C Q0 . 

Replacing E± by E± — £a,±, leaving the first two terms unchanged, we can iterate 
away the last term exactly as above to obtain 

PE ± G z s± W +1 (l - r 2 y + ±W+1 C°° + x »±W+^(l - r 2 Y + ±(X}+1 log(l - r 2 ) + C°°. 
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Repeating this argument proves this proposition. Note that we obtain arbitrarily 
large powers of logarithms, but these correspond to increasingly less singular terms 
in terms of the power s in (1 — r 2 ) s + . □ 

As we would like our operator E± to be localized in the interior of light cone 
(for hyperbolic propagation would spread singularities outside otherwise and E± 
could not satisfy ), it is convenient to consider P a as an operator on tempered 
distributions in 

W£ = {Y: \Y\<1}, 

here equipped with the smooth structure arising from adjoining -Jl — \Y\ 2 to the 
smooth structure induced from the front face (this is what the subscript 1/2 de- 
notes). Let v = (1 — r 2 ) 1 / 2 be a defining function for dM"J^ ■ If only even powers 
of v occur as coefficients of products of vD u and vV, V a vector field on cffi"^ 1 
extended to a neighborhood using the polar coordinate decompositions, then one 
calls the corresponding differential operator even, see [5]. Note that the subspace 
of even elements of C°°(B^ 1 ) is exactly C^iW 1 ^ 1 ). Then: 

Lemma 7.9. P a 6 v~ 2 Diffj^B™^ 1 ) is elliptic and even. 

For a real with A + a 2 — a(n — 1) > 0, —P a is positive with repect to the 
i 2 (B"- x , (1 - v 2 )( n -VI 2 v 1+2 ' J dvduj) inner product on 

vH%(B^, (1 - v 2 )^l 2 v x+2 ° dvduj), 
with duj denoting the standard measure on the unit sphere. 

Proof. As P a is a differential operator with smooth coefficients on all of ff , elliptic 
for r < 1, we only need to analyze its behavior near r = 1. For this purpose it 
is convenient to use the boundary defining function v on B™^ 1 . A straightforward 

calculation using (1 — r 2 ) 1 / 2 !?,. = —(1 — ^ 2 ) 1 / 2 _D !y gives that in fact 

-P a = (D v + i(2a - l)^ 1 + »(n - 3)i/(l - f 2 ) _1 )(l - v 2 )D v 

+ - J A^ + X + a 2 - a(n - 1) 

1 — vr 

-if, , , Un — 3 V 2 , , o, , 

= v~ ( {vD v + i(2o - 1) + \_J 2 )(1 ~ v ){vD v - i) 

v 2 \ 
+y—^2 A w + (A + o- 2 - a(n - l))is 2 J v~ x 

from which the first claim follows immediately. For the second claim we merely 
need to notice that the formal adjoint of D v v = vD v — i with respect to fv^ 1 dv dtu 
f = (l_z,2)(n-3)/ 2l ,2+2 ^ ig /-!(„£>,, _»)/ = vD v +i{2a - 1) + i(n - i)v 2 (1 - v 2 Y x , 
so 

(ti, - 11(1 - ^^AHI 2 + 11(1 - ^ 2 )- 1/2 rf^|| 2 + (A + a 2 - a(n - l))\\u\\ 2 . 

□ 

In fact, it is also convenient to identify the interior of B™^ 1 with the Poincare 
ball model of hyperbolic (n — l)-space H™ -1 using polar coordinates around the 
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origin, letting coshp — v , p is the distance from the origin. The Laplacian on 
jjn-i m these coordinates is 

A H »-i = Dp- i(n - 2) cothpL> p + (sinhp) _2 A w . 

Lemma 7.10. Let s be such that 2s = a — J. T/ien 

(1 - r 2 r s P CT (l - r 2 ) s = vV-'Prif-* 

-v- 1 I Aion-i + cr 2 



= — coshp I Agn-i + a 2 — ^— - — J +(coshp) 4 — -J J cosh p. 

Thus, this conjugate of P a is essentially a compact perturbation of the hyperbolic 
Laplacian, shifted by the eigenparameter (n — 2) 2 /4 — cr 2 . Note that the spectrum 
of Ajjn-i on L 2 (H™ _1 ) is \(n — 2) 2 /4, oo). In fact, we have the following result of 
Mazzeo and Melrose [S]: 

Lemma 7.11. The operator 

2 




L a = A H „-i + a 2 - y—r^ ! + f - I A - 

is invertible on 

L 2 (W n -\p Mn -i) = L 2 (B™^ 1 , (1 - I/ 2j(n-3)/2 I/ i-n d ^ = L a (E^, (sinh p)"- 1 dp) 

/or cr 2 ^ R, and it is Fredholm in a 2 for a 2 G C \ [0, oo). 

For cr > 0, any element of the L 2 -nullspace of L a lies in z/(" _2) / 2+cr C°°(B" / ^ 1 ). 

The inverse L~ x is meromorphic for cr 2 G C \ [0, oo) with finite rank residues, 
maps u k HQ l (M"j 2 , /i H n-i) — > v k H™ + (M™7 2 , /in™- 1 ) continuously, provided that 
\k\ < Recr|. For k > | Recr|, it maps 

v k H™ (B™^, 1 , /%«- 1 ) -> ^^^(B^^Hn-i) + ^("- 2 ^ 2+ff C 00 (B^ 1 ). 

In fact, L^ 1 , defined at first in Recr > 0, extends meromorphically to all of C 
(i.e. the Riemann surface of cr 2 ), as shown in [9j with improvements in [5]: 

Lemma 7.12. The operator L^ 1 defined at first for Recr > as the inverse of L a , 
extends to a meromorphic family of operators 

Ro(a) : ^^(B^S/iHn-O -► ^• J ff"+ 2 (B^ 2 1 , AiH „- 1 ) + ^ ( "" 2)/2+ff C°°(B"- 1 ), 
k > | Recr| with no poles for a ^ pwre imaginary, which satisfies L a Ro(a) = Id 

Moreover, a is a pole of Ro, then L a u = has a non-zero solution 

u G zy ( "" 2)/2+,T C 0O (B n " 1 ). 

Corollary 7.13. For a 2 G C \ [0, oo), Recr > 0, the operator P a is Fredholm, of 
index 0, as a map 

P<r : ^^(^Smh-O - ^-X ro+2 X»i7 2 \MH.»-0 

/or — < k < 2 Recr — ^r', P& 1 is meromorphic, with finite rank poles, and all 
poles satisfy a 2 G K. 
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Moreover, for a > 0, elements of the nullspace of P a on v k H™(W^ 2 , /Xh™- 1 )? h 
1 > 

/2 



as above, lie in ^ 2ct C°°(B i 1 1 7 2 1 ) 



In addition, for k > 2Re<7 — ? whenever P a is invertible on L 2 , 

P- 1 : v k H™{WlJ 2 \^-,) - v k - 2 H™ + *{WlJ 2 \n m -,) +^C°°(B"- 1 ). 
Finally, R(cr) — P~ , Re a > ? extends to a meromorphic family 

jn-l \ , 2 rrm+2 /-iron— 1 \ , 2a noc /"iron— 1 ^ 



A: > 2| Rca| — lyzift no po/es /or cr 7^ pure imaginary, and P a R(a) = Id on 

nn— 1 
6 l/2 



i/kjjm^n . Hffln-i^ k as above. If a is a pole of Rq, then P a u = has a non-zero 



solution 

Proof. P a = -v cr -9- 1 L <T v-' T+ 9- 1 , so 

p-l _ +i i -i I/ - -+f +i_ 

□ 

Note that 1 just barely fails to be in iH"- 2 )/ 2 !, 2 ^^ 1 , ft^-i), while 

,2tr/-JOO /"iron— 1 \ ^— , ,fc r 2 /iron — 1 



n-2 



= v- x (yD v - - v 2 ){vD v - i)v~ l + -A u + X + a 2 + a{n- 1), 



for k < 2 Recr > . 

If A < 0, s+(A) > 0, and P$ + (X) fails to be invertible on the spaces listed above 
as -P,s + (A)i /2s+< - A ' ) = 0, and iy 2s +( A ) lies in these spaces. However, we claim that P a 
is invertible for a > s + ( A). In fact, 

-P a = -v 2a (y- 2a P a v 1 < y )v-' la = -v 2a P- a v- 2a 

1 

1 ~v 2 

with adjoint taken relative to (1 — v 2 )( n ~^/ 2 v 1 ~ 2,y dv dtu. The first two terms 
are positive with respect to the corresponding I? space, while the roots of A + 
a 2 + a(n — 1) are exactly s±(A), so A + a 2 + a(n — 1) > for a > s+(A). As 
is 2a C 00 {M n - 1 ) C i^B"- 1 ,!/ 1-2 " dvdto), it follows from Corollary \TM that P a has 
no nullspace in the listed spaces, so it is invertible. (A different way of arguing 
would have been to note that vPs + (\)V has a positive eigenfunction, j/~ 1+2s +( A ), 
which thus must correspond to the bottom of the spectrum.) 

That for A > the poles do not occur follows from the following lemma as 

^ 2 (B^ /2 1 ,^- 2Ro > H „- 1 ) = i, 1 -f+ RcCT L 2 (B' 1 7 2 1 , / i H „- 1 ), 

and 1 - f + Re a < 2 Re a - 2=2 if R e a > 0. 

Lemma 7.14. P a satisfies 



Pa = -(A- - q - n(D r - - -) - r~ 2 A u 



^ gx 1 - r 2 ' v ' v 1 - r 2 ' ~ 1 - r 2 

= -(£>„ + iCTi/" 1 )*^ - v 2 )(D v + lav' 1 ) - — A w - a 2 ^ 2 - A 

1 — v 1 

with the (formal) adjoint taken with respect to the measure 

(i = r n - 2 (l - r 2 )- Rca drduj = (1 - V 2 )^v 1 - 2Kca dv duj = V 71 - 2 ™*" fi U n-x. 
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Corollary 7.15. Suppose that A < (n — l) 2 /4. Then P a is invertible for a > 
max(0, s+(A)). 

In fact, we can analyze the poles of the analytic continuation R{o~) rather accu- 
rately using special algebraic properties of P a . Unlike the proceeding considerations, 
which were rather general, i.e. hold for operators of the same form, the following 
relies on the precise form of P a . 

Lemma 7.16. The following identities hold: 

P ff _ 2 Ay = A Y P a , P a + 2V 2 °^A Y v- 2 ° = v 2 ° +i A Y v- 2a P a . 

Proof. First, as Ay is homogeneous of degree —2 with respect to dilations on Y, 
[yd y , A Y ] = -2Ay, so [YD Yl A Y ] = 2iA Y . As 

P a = (YD Y -i(n-l- a)){YD Y + ia) - A Y - A, 

we deduce that 

A Y P a = P a A Y + [Ay, (YD Y f + i{2a - (n - l))YD Y ] 

= P a A Y - 2iA Y (YD Y ) - 2i(YD Y )A Y + 2(2a - (n - l))Ay 

= P a A Y - 4A y - Ai{YD Y )A Y + 2(2a - (n - l))Ay 

= ((YD Y ) 2 + i(2a - (n - l))YD Y + a(n — 1 — a) — Ay - A 

-4 - Ai{YD Y ) + 4cr - 2(n - 1)) Ay 
= ((YD Y ) 2 + i(2(a - 2) - (n - l))YD Y + (a - 2){n -I -a + 2) 

-Ay - A) Ay 

= P*- 2 A Y . 

Thus, using P_ CT = v~ 2 ° ' P a v 2a with a replaced by a + 2 first, then with a replaced 
by -a, 

P a+2 v 2a+i A Y v- 2 ° = v^P.^Ayv- 2 * 

= v 2 °^A Y P_ a v- 2 ° = v 2 ° +4 A Y is- 2a P„ 

as claimed. □ 

Lemma 7.17. Suppose that a is such that P a+2 w = 0, w e v 2 ^ +2 ^> C 00 {W 1 ' 1 ) 
implies w = 0. If P a u = for some u G v 2 " C°° {B> n ~ l ) then either a e s±(A) - N 
or u = 0. 

Proof. If P a u = 0, then by the previous lemma, P CT+2 ^ 2cr+4 Ay^ _2cr ?i = 0. More- 
over, v- 2a u e C°°(B™- 1 ), so w = v 2a+i A Y is- 2a u e i/ 2 ( CT+2 )C°°(B"- 1 ), hence w = 0. 
Thus, v = v~ 2rT u e C°°(B™- 1 ) satisfies Ayt) = and P_ CT t; = i/- 2cr P_ CT i/ 2ff ?; = 
v- 2a P a u = 0. Thus, (P_ CT + A Y )v = 0, so 

((yPy) 2 - i(n - 1 + 2cr)rPy - (A + cr(n - 1 + a))) v = 0. 

Factoring the operator as (YD Y + ia + )(YD Y + ictJ) with 

a± = -^—^ -<r± \j (^y^) ~ x = M A ) - CT > 

we deduce that v satisfies either (YD Y + ia + )v = or (YD Y + ia-)v = 0, i.e. v 
is homogeneous of degree a+ or degree a-. But v is C°° at the origin, so, unless 
v = 0, in either case the corresponding a must be a non-negative integer, i.e. 
s±(A) — a = m E N, so a 6 s±(A) — N, proving the lemma. □ 
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Corollary 7.18. The only possible poles of R(o~) are a £ s±(A) — N. In particular, 
if m is a positive integer, R(cr) is regular at a = s±(A) +m unless s+(A) — s_(A) = 

Proof. As noted in Corollary |7.131 a is a pole of R if and only if there exists a non- 
zero u € i/ 2ct C 00 (B™" 1 ) such that P a u = 0. Moreover, if Re a > C, C sufficiently 
large (depending on A) , then there exist no such non-trivial u by Corollary 17.151 
Correspondingly, if Re a G (C — 2, C] and a is a pole of R, then the previous lemma 
shows that a G s±(X) — N. Proceeding inductively we deduce the corollary. □ 

Now, if a is not a pole of R, then given / G C°°(B™^ 2 1 ), P a v — f can be solved 
with v G iy 2a C oc (M n - 1 ). 

If Reer > and we extend v as to the rest of the fiber of the front face over y', 
P a v is thus the extension of /. 

In fact, as long as 2a ^ — N+, we can extend v by expanding in Taylor series to 
finite order, v — J2j=o p2 ' ' a j + v 2N+2 v' ', v' C°° near 9B" -1 . If we choose N large 
enough so that 2Recr 4- 2N + 2 > 0, we can extend v 2a v' to ff by extending it as 
0. On the other hand, we can extend v 2a+2 ^ aj as (1 — |y| 2 )^. Thus, we obtain 
a distribution v on ff. Now P a is a second order differential operator with C°° 
coefficients, so P a (l - \Y\ 2 f + +J has the form (1 - \Y\ 2 ) a +^ 2 b' j , with b' 3 smooth, 
and as the principal symbol of P a vanishes on the conormal bundle of it 
in fact has the form (1 — \Y\ 2 ) a t + -'~ 1 bj, with bj smooth, as long as a + j is not a 
non-positive integer. In particular, we deduce that P a v = provided that P a v = 0. 

This is the argument that requires using the analytic extension of R to Re a < 0, 
which gives solutions v £ i/ 2(7 C° (M n ^ 1 ) rather then using solutions involving the 
other indicial root, 0, which would give rise to v 6 C°°(B n_1 ), and hence allow P a v 
to have delta distribution terms at 9B™ _1 . In particular, for Reer < 0, we cannot 
simply use the conjugate (in the sense of Lemma 17. 10p of L_„. 

If a e s_(A) + N+, 2a e -N+ can hold only if 2y/(n - l) 2 /4 - A G N+; it can 
never hold if a G s+(A) + N+. We thus deduce that with s = s+(A), or s = s_(A) 
under the additional assumption that s+(A) — s-(A) ^ N, we can solve away the 
error in Taylor series to obtain 

2n -\- 1 

E ± G x s I m ^([X x y + ;diag y+ ],n m(s) = -s - — s = J±(A), 

with E ± - E .± G x s+1 P n ^([X x Y + ;<tia,g Y+ ],F), and E± supported inside the 

light cone, PE G C°°(X x Y+). This remaining error can be removed using the 
results of Section [3] to obtain the same conclusion with PE± — near Y + . The 
standard FIO contruction allows one to obtain E± with the same properties, except 
PE± supported near Y_, vanishing in a neighborhood of K_. We have thus proved: 

Proposition 7.19. Suppose that s — s+(A), or s — s_(A) under the additional 
assumption that s+(A) — s_(A) ^ N, i.e. A ^ ^ 4 ~ m , m G N. Then there exists 

E ± (X) G x s T m ^([X x y+;diag y+ ],F), m(a) - -s - ^±1 

satisfying PE±(\) = near y + x y + , -E±(A) supported inside the light cone near 
Y + x y + . anrf 

E±(X)<t) = x s±w v, veC°°(X), v\ Y+ =4> 



THE WAVE EQUATION ON ASYMPTOTICALLY DE SITTER-LIKE SPACES 



37 



for all <f> £ C°°(Y+). Moreover, a m ( s )(E±) never vanishes. 

We let 25(A) = E + (X) © 25_ (A) be the Poisson operator (near Y + , where it solves 
PE±(X) = 0). However, much as it is useful in the interior of X to renormalize using 
powers of the Laplacian, the same holds here. The renormalization depends on the 
choice of x modulo x 2 C°°(X). So let Aj, denote the Laplacian of the boundary 
metric h, define A' h analogously to the case of Cauchy surfaces, i.e. is Id on the 
nullspace of A/j, and is A/, on its orthocomplement. The renormalized Poisson 
operator is then 

25(A) =^ + (A)(A', l )( s +( A )-™/ 2 )/ 2 ©25_(A)(A' fc )( s -W-™/ 2 ^ 2 . 

The n/2 in the exponent of A^ is somewhat arbitrary, it is used to normalize 
FIO's below to be zeroth order; any quantity differing from s±(A) by a constant 
(s-independent) amount would work. By Proposition 17.191 the two components 
of 25(A) lie in a; 1 -"/ 2 2- 5 / 4 ([X x Y+; diag y+ ], F; (C 2 )*), i.e. they have the same 
regularity in the interior of X x Y + as well as the same behavior at the boundary. 

Proposition 7.20. Suppose s+(A) - s_(A) N, i.e. A ^ (n-if-m 2 ^ m g N For 
to sufficiently close to I, the map sending scattering data at Y + to Cauchy data at 
St given by 

S+j : C°°(Y + ) 2 3 [g + ,gj) » (u\ Sto ,d x u\sJ € C°°(S to ) 2 , 

where u is the solution of Pu = given by Theorem 15.51 is the Fourier integral 
operator with Schwartz kernel 2?(A)|s + ( e ) x y + © 9 a 2?(A)|E + (e)xY+- 
Moreover, the renormalized map 

S +M = R t0 E(X) =R t0 E + (\)(A'J s +^- n /V/ 2 (B R t0 E-(\)(A' h Y s -W- n ^/ 2 
e I°(S t0 x Y+, F n (S t0 x Y+); £(C 2 , C 2 )), 

with Rt being the Cauchy data map at to, u i— > ((A( o ) 1 / 2 it|5 t ,Vu\s t ),V a vector 
field transversal to St , is an invertible Fourier integral operator. 

Proof. Let t\ < to, but still sufficiently close to 1. Let x G C°°(X) be identically 1 
in a neighborhood of T > to, supported in T > t\. Let u be the solution of Pu = 
given by Theorem 15.51 and let v — J Y \E+{X)g+ dy + J Y \E-(\)g- dy. Then at 
Y + , v has the asymptotics required by Theorem 15. 51 and Pv = [P, x]{Jy E+g+ dy + 
f y E_g_dy) is supported where T £ (ti,t ). For I > max(i,/(A)), let v\ £ 
Hq' 1 '~ 1 (X) be the solution of Pvi — —Pv, As in the proof of Theorem 16. 1[ v% 
is identically for T > to, is in C°°(X°), and has an asymptotic expansion at F_ as 
in Theorem 15. 51 Thus, v + v± has all the properties of u required by the uniqueness 
part of Theorem 15.51 so 

u = J xE+g+ d V + J^ X E -g- dy + «i , 
and v\ = at St . Thus, *S+,t indeed has Schwartz kernel 

E W\s tQ xv+ ffi d x E(X)\ St0 xv+, 

and we have a similar expression for E(X). As F is transversal to St and the 
restriction map to S to is a Fourier integral operator of order 1/4, 25(A) |,s to x y+ is an 
FIO of order —1, while R to E(X)\s to xy + is an FIO of order 0. 
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In order to prove the invertibility of 5+,t , it suffices to show that it is elliptic in 
the sense that tQ S+.t a and <S+,t «S+ t are elliptic pseudo-differential operators, 
where the adjoint is taken with respect to the Ricmannian densities on St and Y + . 
Once this is shown, it follows that both the nullspace of iS+,t and of its adjoint 
must lie in smooth matrix-valued functions, and are finite dimensional. Consider 
for instance <S+,t - F° r such smooth Cauchy data (17+, gJ) at Y + , the corresponding 
solution of Du = is smooth in X°, of the form given by Theorem 15. 5[ and the 
vanishing of its Cauchy data at St implies that in fact u vanishes identically, hence 
g± = 0, so <S+,t has trivial nullspace. On the other hand, suppose that to is 
not injective, i.e. S+j a is not surjective (e.g. on the L 2 -spaces). Any element of the 
nullspace of <S+ tg is smooth, so in this case there exist smooth non-zero Cauchy 
data (ifjQjtpi) & t St which are L 2 -orthogonal to the range of 5+,t . Let u be the 
solution of Pu — with these Cauchy data. Let (g+ ,g~) be the leading coefficients 
of the asymptotics at 1+, as in Theorem 15. 51 Then u = E+(\)g+ +£'_(A)g_ (since 
the right hand side has the same asymptotics at Y+ as the left hand side, so they are 
equal by the uniqueness part of Theorem [53]) . Therefore (i/>o>'0i) are in the range 
of <S+,t , so they vanish, which gives a contradiction. Thus, t is also injective. 
This proves the invertibility of <S+,t given its ellipticity. 

In order to prove ellipticity, one needs to compute the principal symbol of 
. t <S+,t and S +tto S^_ t . Consider first the latter. For each a = (z, C) 6 T*S to 
there are two bicharacteristics of □ which contain a point over z S St whose im- 
age in T*S ta = Tg t ^X/N*St is (z,(). Let the corresponding points in T*S to 
be (Xj — {to,Z,£j,Oi i = ; where £ is the dual variable of the first co- 
ordinate, T. These bicharacteristics emanate from S±Y + (one from S1Y+, one 
from S*lY + ): let (3j = (yj,r)j), j = +,— , be the corresponding points. Let 
E± = E±(A' h )( 8± W~ n / 2 " 2 . Then Hormander's theorem on the composition of 
FIO's shows that the principal symbol of S+ t t Sl. t at a = (z,C) is a constant 
times 



E 



^{a)a{E + )( aj ,P J ) aiA'^iaME^a^fr) 



a(V)(a 3 )a(E-)( aj ,P 3 ) 



E( 



a(A' h )y 2 (a) a{E ± )(a il 0j) < 7{V){a j ) a(E + )(a j ,0 j ) 
a(E + )\ 2 + \a(E^)a(A') 



a(V')(a J -)(7(/v )(a ; -,a ; ). 



where 



_ a{V) 



, and where on the right hand side the various principal sym- 



bols are evaluated at the same points as on the left hand side, but suppressed in 

r 1 r ■ 1 

notation. Thus, the principal symbol has the form J^. -Cj _ , J , 2 with Cj > 0, 

[ r j \ r J\ \ 

and a straightforward calculation shows that this matrix is positive definite, hence 
invertible, provided r + ^ r_. But r + = r_ would imply that a + = a_, which is 
not the case, so we conclude that S-\- t t S^_ to is indeed elliptic. 

The calculation for 5+ t «ST t is similar. In this case, for each j3 = (y, rf) G S*Y + , 
there are two corresponding bicharacteristics, again one including a point in S1Y + 
and one in S*L Y+, which then cross Tg t X at a.j — (xj, Zj, S,j,Tjj), j — +, — ■ Thus, 
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the principal symbol at /3 = {y, rj) is 



E 

3 



a{U h )V*{a j )a{E + ){ aj ,(3) a(V) («>(£+) (ay , (3) 



a(A' ft )V2 (aj)a(E- ) (ay , p) a(V) (ay )a(E-) (ay , (3) 



a(V)( aj )a(E + )( aj ,P) 
J2(v(&') + W(V)\ 2 )W(E + )\ 2 



a(V){ aj )a(E-){a 3 ,(3) 



where now ry = ^j^ri > and where again on the right hand side the various principal 
symbols are evaluated at the same points as on the left hand side, but suppressed 
in notation. 

Now a_su(E+) and cr_K,u(E-) satisfy the same first order linear ODE along 
bicharacteristics, so their ratio along each bicharacteristic is constant, hence are 
equal to the ratio evaluated at the 'initial point' at the front face of [X x Y+; diagy + ] 

(where cr(E + ) has to be replaced by a(x n ^ 2 ^ 1 E + ), etc.). For a given (y',r)'), the 
projection of the two bicharacteristics hit the front face at (y',Y), Y = ±77', and 
the bicharacteristics themselves hit the cotangent bundle over the front face inside 



N{ y , X) F at 



-yd|y|. Thus, 



C r(x"/ 2 - 1 £_)(y',y il -Y J -d|y| 



aixo/t-iE+W^-YjdW) 

But these can be calculated from the normal operators, which are explicit, hence 
are easily evaluated as 1, resp. e t7: ( s +( xs, ~ s -( x ^ . Thus, if s+(A) — s_(A) is not an 
even integer, which we are assuming, the ry are unequal, so 5^ t <S+,t is indeed 
elliptic, finishing the proof. □ 

We are now ready to prove one of our main results, that the scattering operator 
is a Fourier integral operator. 

(n-l) 2 -m 2 



Theorem 7.21. Suppose that s+(A) - s_(A) ^ N, i.e. A ^ 



€ N. 



Then 5(A) is a Fourier integral operator with canonical relation given by S c \, and 
S = 5(A) is an invertible elliptic Oth order Fourier integral operator with the same 
canonical relation. 



Proof. This is immediate from S = S_ 



-l+e 



o Ci_ c -i+e ° S+ i-e for e > small. 



Indeed, all operators are Fourier integral operators by Proposition 17.201 (applied 
also at YL) and Proposition 16.21 with canonical relation given by the appropriate 
restriction of the bicharacteristic flow. Thus, the projection of the canonical relation 
to each factor for each of them has surjective differential, so the composition is 
transversal, and Hormander's theorem can be applied. As 

5(A) 

= ((A^)" s +W/ 2 +"/ 4 © (A' h )- s -^ +n / i )S(\)((A' h ) s +^/ 2 - n / i © (A'^-W/ 2 -"/ 4 ), 
and the first and last operators are pseudodifferential, the theorem follows. □ 



Theorem 11.41 follows similarly, as the propagator mapping Cauchy data at dif- 
ferent T-slices to each other is an invertible FIO, so it suffices to consider the case 
to close to 1, in which case the inverse given by Proposition 17.201 proves the claim. 
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